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ABSTRACT 


In  urbanized  terrain,  radiowave  propagation  is  subjected  to  fading  on  large-scales 
and  small-scales  that  would  impede  on  the  quality  and  reliability  of  data  link  transmis¬ 
sion.  This  would  have  implications  in  many  military  applications.  One  example  is  the 
performance  of  unmanned  aerial  vehicle  (UAV)  data  and  communications  links  in  com¬ 
plex  urban  environments. 

The  purpose  of  this  research  is  to  study  the  effectiveness  of  diversity  techniques 
on  the  performance  of  urban  UAV  data  and  communications  links.  The  techniques  inves¬ 
tigated  were  spatial,  polarization,  and  angle  diversities. 

The  ray  tracing  software,  Urbana  Wireless  Toolset,  was  used  in  the  modeling  and 
simulation  process.  The  various  combinations  of  diversity  techniques  were  simulated  us¬ 
ing  a  realistic  urban  city  model.  For  the  few  transmit-receive  geometries  examined,  it  was 
found  that  angle  diversity  with  a  directive  antenna  provided  the  greatest  increase  in  signal 
strength  relative  to  the  no  diversity  case. 
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EXECUTIVE  SUMMARY 


In  military  operations  in  urbanized  terrain  (MOUT),  achieving  reliable  tactical 
communications  can  be  challenging.  In  urban  areas,  the  buildings  and  structures  can  re¬ 
sult  in  multipath  that  leads  to  fading  and  excess  path  loss  in  radiowave  propagation.  This 
problem  is  evident  especially  for  mobile  land  forces  (i.e.,  troops  and  vehicles)  that  are 
dispersed  at  street  level  confined  on  two  sides  by  high-rise  buildings. 

An  unmanned  aerial  vehicle  (UAV)  with  a  communications  payload  can  be  con¬ 
figured  as  airborne  communications  node  (ACN)  for  urban  communications  capability.  A 
UAV  that  operates  as  a  communication  relay  for  mobile  units  operating  deep  in  urban 
canyons  can  shorten  the  communications  distance,  thus  reducing  path  loss  and  overcome 
non-line-of-sight  (NLOS)  situations.  However,  multipath  and  signal  fading  are  still  sig¬ 
nificant. 

The  ray  tracing  software,  Urbana  Wireless  Toolkit,  from  Science  Applications  In¬ 
ternational  Corporation  (SAIC),  was  used  to  model  and  simulate  the  effects  of  diversity 
on  urban  radiowave  links.  It  is  a  computational  electromagnetic  (CEM)  tool  utilizing  a 
ray  tracing  engine  that  combines  physical  optics,  geometric  optics,  and  diffraction  phys¬ 
ics  to  produce  a  3-D  simulation  for  radiowave  propagation. 

For  urban  radiowave  propagation  modeling,  the  input  parameters  for  Urbana  are 
CAD  facet  models  for  the  buildings  and  roads,  including  its  surface  material  (e.g.  con¬ 
crete,  earth,  glass  and  dielectric  properties).  The  other  inputs  are  the  placement,  strength, 
and  antenna  patterns  of  transmitters  and  receivers.  The  outputs  are  composite  signal 
strength,  angle  of  arrival  and  time  delay  of  each  arriving  signal  ray  at  each  receiver  or  ob¬ 
servation  point. 

In  this  thesis,  three  diversity  techniques  were  investigated  to  assess  their  perform¬ 
ance  in  an  urban  environment.  They  were  (1)  spatial  diversity,  (2)  polarization  diversity, 
and  (3)  angle  diversity.  These  techniques  were  simulated  for  a  UAV  platform  operating 
over  a  small  city  model.  The  findings  can  be  summarized  as  follows.  Angle  diversity  has 
shown  the  most  significant  improvement  in  received  signal  strength  as  compared  to  the 

xvii 


other  two  techniques,  but  the  implementation  of  such  directional  antennas  would  be  com¬ 
plex.  Spatial  diversity  with  two  antennas  on  receive  has  shown  better  overall  improve¬ 
ment  as  compared  to  polarization  diversity.  Further  improvement  is  possible  with  more 
spatially  uncorrelated  antennas,  but  its  implementation  could  be  limited  by  the  availabil¬ 
ity  of  space.  Polarization  diversity  on  transmit  and  receive  has  shown  good  received  sig¬ 
nal  strength  only  in  a  few  cases  over  the  spatial  diversity  technique.  However,  a  combina¬ 
tion  of  polarization  diversity  on  transmit  and  spatial  diversity  on  receive  has  better  over¬ 
all  perfonnance.  An  advantage  of  polarization  diversity  is  that  the  two  orthogonal  linear 
polarizations  can  usually  be  obtained  with  a  single  antenna.  This  is  especially  appealing 
for  UAV  applications,  where  space  on  the  platform  is  limited. 

This  research  has  demonstrated  the  potential  of  implementing  diversity  techniques 
for  improving  UAV  data  links  in  MOUT.  It  should  be  noted  that  the  simulation  per¬ 
formed  did  not  cover  all  possible  system  and  environment  parameters.  A  limited  number 
of  transmitter-receiver  geometries  were  examined.  However,  care  was  taken  to  select 
cases  that  included  direct  line-of-sight  (strong  signals)  and  non-line-of-sight  in  shadows 
(faded  signals).  This  allows  a  comparison  of  various  techniques  under  a  range  of  condi¬ 
tions  so  that  common  general  conclusions  could  be  reached. 


I.  INTRODUCTION 


A.  UAV  IN  MOUT 

Military  operations  in  an  urban  environment  is  a  tricky  business,  as  history  has 
shown.  From  Stalingrad  (1944)  to  Somalia  (1991)  to  the  current  Operation  Iraqi  Free¬ 
dom,  there  are  many  examples  of  complex  military  operations  in  urbanized  terrain 
(MOUT)  where  conventional  battlefield  tactics  do  not  apply.  Where  urban  operations 
cannot  be  avoided,  some  authorities  [1]  have  advocated  aerospace  operations  to  be  the 
key  for  successful  MOUT.  One  increasingly  important  aspect  of  aerospace  operations  is 
the  application  of  unmanned  aerial  vehicles  (UAV).  Figure  1  shows  some  uses  of  UAVs 
for  intelligence,  surveillance,  and  reconnaissance  operations,  and  as  well  for  communica¬ 
tions  relay  node. 


RAMD  W?)l &J-5 


NOTE:  BM  =  balllespace  management;  C2  =  command  and  control;  COM  = 
communications;  I IA I !  =  High  Altitude  Endurance;  IMINT  =  image  intelligence; 
SIGINT  =  signals  intelligence. 


Figure  1.  UAV  sensor  communications  concept  (From  Ref.  [1].) 


The  suitability  of  UAVs  for  MOUT  is  mainly  due  to  the  fact  that  they  can  be 
small  and  deployed  near  potential  targets  with  a  low  probability  of  detection.  They  can 
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have  long  endurance,  which  is  important  for  surveillance  missions  that  could  last  for  days 

[2]. 


In  recent  years,  there  is  an  emergence  of  mini-UAVs  for  use  in  MOUT.  Figure  2 
shows  a  vertical  take-off  and  landing  (VTOL)  Cypher  UAV  in  operation  in  a  confined 
area.  The  idea  of  a  small  UAV  that  can  be  used  to  provide  overhead  surveillance,  remote 
sensing,  and  communications  relay  would  be  very  useful  to  troops  on  the  ground. 


Figure  2.  VTOL  in  MOUT  (From  Ref.  [3].) 

A  mini-UAV  can  enhance  situational  awareness  by  effectively  extending  the 
troop’s  visual  range  by  allowing  views  around  buildings,  terrain  and  other  obstacles  as 
shown  in  Figure  3.  Thus,  reconnaissance,  surveillance  and  target  acquisition  can  be  per¬ 
formed  without  exposing  troops  to  danger  [3].  Similarly,  communications  effectiveness 
can  be  enhanced  by  using  mini-UAVs  as  relay  nodes  to  reach  out  to  troops  in  non-line- 
of-sight  (NLOS)  situations. 


Figure  3.  VTOL  providing  situational  awareness  (From  Ref.  [3].) 
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Table  1  shows  some  examples  of  production  mini-UAVs  currently  on  the  market. 
The  possible  sensors  and  data  link  payloads  are  listed.  The  data  link  specifications  pro¬ 
vide  useful  information  that  were  adopted  in  the  modeling  and  simulation  performed  dur¬ 
ing  the  course  of  this  research. 


NAME  OF 
UAV 

TYPE 

MANUFACTURER 

SENSORS 

DATA  LINK 

Vigilante  496 

Rotary 

SAIC 

EO/IR,  SAR, 

ESM 

900MHz, 

2.4GHz 

Phoenix 

Wing 

BAE  Systems 

IR 

J-band 

Hermes  450S 

Wing 

Elbit  Systems, 

Silver  Arrow 

EO,  SAR, 

Comms  relay 

C/L-band 

Dragon 

Drone 

Wing 

BAI  Aerosystems 

Video/TV 

UHF  uplink,  L- 
band  video  with 
data  sideband 
downlink 

Heliot 

Rotary 

CAC  Systemes, 

Dragon  Fly,  EDT 

CCD,  IR 

S/Ku-band 

Table  1.  Examples  of  mini-UAVs  (Compiled  from  Ref.  [4,  5].) 


B.  IMPROVING  TACTICAL  COMMUNICATIONS  IN  MOUT 

When  UAV  sensor  payloads  are  used  to  gather  important  intelligence  and  surveil¬ 
lance  information,  the  ability  of  the  UAV  to  communicate  reliably  via  wireless  data  links 
to  ground  units  is  crucial.  In  an  urban  environment,  radio  signals  can  be  severely  de¬ 
graded  due  to  multipath  and  fading,  and  this  can  seriously  impact  the  performance  of  a 
UAV  data  link. 

Mobile  land  forces  (i.e.,  troops  and  vehicles)  in  MOUT  are  usually  dispersed  in 
so-called  urban  canyons;  that  is,  at  street  level  confined  on  two  sides  by  high-rise  build¬ 
ings.  Each  dispersed  unit  would  need  to  maintain  communications  links  among  each 
other  and  with  headquarters,  for  sharing  of  information  to  ensure  real-time  situation 
awareness.  These  units  are  mobile  and  can  be  deeply  embedded  in  the  urban  canyons.  As 
such,  it  would  be  difficult  to  maintain  effective  communications  due  to  the  effects  of 
multipath  and  fading,  excess  path  loss,  and  NLOS  situations. 
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In  order  to  alleviate  some  of  the  potential  communications  problems,  besides  pro¬ 
posing  operational  limitations,  a  potential  solution  is  the  use  of  newly  evolving  technolo¬ 
gies  including  software  radios,  ultra-wideband  waveforms,  antenna  diversity,  and  UAV 
relays  [6]. 

A  UAV  with  a  communications  payload  can  be  configured  as  airborne  communi¬ 
cations  node  (ACN)  for  NLOS  communications  capability  between  mobile  land  forces 
and  headquarters  [7].  A  UAV  can  reduce  the  communications  link  distance  and  operate 
as  a  communication  relay  for  mobile  units  operating  deep  in  urban  canyons.  Diversity 
techniques  (frequency,  space,  polarization,  and  angle)  can  be  implemented  on  both  the 
mobile  land  units  and  UAV  to  improve  the  degraded  radio  signals  which  are  subjected  to 
multipath  and  fading. 


C.  OBJECTIVES  OF  THESIS 

The  purpose  of  this  research  is  to  investigate  how  to  exploit  diversity  techniques 
in  mitigating  the  effects  of  multipath  and  fading.  The  focus  is  to  seek  improvement  on  the 
performance  of  UAV  data  links  in  urban  environments.  The  techniques  to  be  explored  are 
spatial,  polarization,  and  angular  diversities. 

This  research  involves  modeling,  simulation,  and  analysis  of  the  UAV  data  link 
performance.  Various  spatial,  polarization  and  angular  diversity  techniques  were  imple¬ 
mented  in  simulations  in  order  to  study  the  effects  on  signal  distribution  in  an  urban  envi¬ 
ronment.  The  modeling  and  simulation  provides  insight  into  the  characteristics  of  urban 
radiowave  propagation  and  the  effectiveness  of  diversity  techniques  on  the  link  perform¬ 
ance.  The  mitigating  effects  on  radio  signal  fading  are  demonstrated.  In  summary,  this 
investigation  suggests  possible  methods  of  improving  the  perfonnance  of  UAV  data  links 
operating  in  urban  environments. 

Parsons  [8]  highlighted  that  propagation  methods  based  on  the  deterministic  ray 
theory  have  enormous  potential.  Ray  techniques  have  been  used  extensively  in  many  in¬ 
vestigations  in  recent  years,  for  both  indoor  and  outdoor  environments.  Similarly,  ray 
tracing  was  the  primary  investigative  method  for  the  modeling  and  simulations  of  the 
problems  in  this  thesis.  The  ray  tracing  software  Urbana  from  SAIC  was  used.  In  addi- 
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tion,  Matlab  software  was  used  to  assist  in  the  post-processing  and  visualization  of  simu¬ 
lation  data. 


D.  THESIS  OUTLINE 

Chapter  II  discusses  the  characteristics  of  the  urban  radiowave  propagation  chan¬ 
nel,  and  presents  some  commonly  used  empirical  models. 

Chapter  III  reviews  various  topics  in  spatial,  polarization  and  angle  diversity  tech¬ 
niques,  and  discusses  how  they  can  be  used  to  improve  data  link  performance  in  an  urban 
environment. 

Chapter  IV  describes  the  modeling  process  using  the  Urbana  software  and  how 
operation  of  a  UAV  data  link  with  diversity  can  be  simulated. 

Chapter  V  describes  the  simulations  and  presents  data  for  selected  cases. 

Chapter  VI  presents  some  conclusions  and  suggestions  for  further  research. 

Appendix  A  lists  the  Urbana  codes  used  in  the  modeling  and  simulations.  Appen¬ 
dix  B  compiles  the  Matlab  codes  used  in  this  thesis.  Appendix  C  presents  the  coordinates 
system  used  in  Urbana,  and  Appendix  D  lists  the  additional  simulations  and  post¬ 
processing  results. 


5 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


6 


II.  CHARACTERIZATION  OF  URBAN  RADIOWAVE 
PROPAGATION  CHANNEL 


Urban  radiowave  propagation  is  not  a  single  propagation  mechanism  but  a  com¬ 
plex  combination  of  direct,  reflected,  refracted,  and  diffracted  wave  components.  De¬ 
pending  on  the  relative  strengths  and  phases  of  these  components,  they  can  interact  de¬ 
structively  leading  to  fading  [9]. 

Propagation  between  a  transmitter  and  receiver  is  subjected  to  three  scales  of  sig¬ 
nal  fading  as  shown  in  Figure  4.  They  are:  (1)  large-scale  path  loss,  which  refers  to  the 
signal  variations  due  to  the  spreading  of  wavefront  in  space  over  the  given  environment; 
(2)  shadow  fading,  which  is  the  signal  variation  due  to  the  terrain  features;  and  (3)  small- 
scale  fading,  which  are  signal  variations  due  to  multipath  between  the  transmitter  or  re¬ 
ceiver  [10].  In  this  chapter,  these  three  types  of  fading  are  discussed. 


Figure  4.  Received  signal  in  a  multipath  environment  (After  Ref.  [10].) 


A.  LARGE-SCALE  PATH  LOSS  MODELS 

There  are  many  empirical  models  that  describe  the  characteristics  of  large-scale 
signal  loss  in  an  urban  environment.  However,  empirical  models  are  area  or  city  specific, 
which  might  not  be  relevant  if  the  city  of  interest  differs  from  the  reference  model.  Be¬ 
sides  the  empirical  models,  there  are  deterministic  physical  models  which  provide  for  the 
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flexibility  of  prediction  for  specific  areas  or  cities.  Ray  tracing  techniques,  which  were 
used  in  this  thesis,  are  an  example.  A  few  of  the  large-scale  path  loss  models  are  re¬ 
viewed  in  the  following  paragraphs. 


1.  The  Okumura-Hata  Model 

This  is  an  empirical  model  based  on  measurements  done  in  Tokyo  for  frequencies 
between  200  MHz  and  2  GHz.  Its  urban  area  is  defined  as  a  built-up  city  or  large  town 
with  large  buildings  and  houses  having  two  or  more  stories.  For  medium  to  small  cities, 


the  total  loss  is  given  as  [11] 

L[dB]  =  A  +  B\ogR-E  (1) 

where 

A  =  69.55  +  26. 16 log  fc  - 13.82  log hb,  (2) 

B  =  44.9  -6.55  log (3) 
£  -(1.1  log /c  -0.7)/im  -(1.56log/c  -0.8),  (4) 


and  Equation  (4)  applies  to  medium  to  small  cities. 

This  model  is  valid  for  carrier  frequencies  in  the  range  of  150  MHz  <  fc  <  1.5 
GHz,  transmitter  antenna  heights  in  the  range  of  30  m  <  hh  <  200  m,  receiver  antenna 
heights  in  the  range  of  1  m  1 0  m,  and  a  transmitter  receiver  spacing,  R  >  1  km. 

Under  this  model,  there  are  other  equations  for  large  cities,  suburban  and  open  areas, 
which  are  not  presented  here  since  the  city  model  of  interest  in  this  thesis  falls  into  the 
category  of  a  small  city. 

2.  The  COST  231-Hata  Model 

In  this  model,  the  Okumura-Hata  model  for  medium  to  small  cities  is  extended  to 
the  frequency  band  1.5  GHz  <  fc  <  2  GHz.  The  total  loss  is  given  as  [11] 

T[dB]  =  F  +  51ogR-£  +  G  (5) 
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where 


F  =  46.3 +  33.9 log /c  -13.82 log and 

fO  dB,  medium-sized  cities  and  suburban  areas 
G  =  \ 

[3  dB,  metropolitan  areas. 


(6) 

(V) 


3.  Free  Space  Loss 

In  the  physical  models  that  are  described  next,  the  loss  is  given  relative  to  the  free 
space  path  loss,  and  is  referred  to  as  excess  loss.  The  equation  for  the  free  space  loss  is 
given  as  [1 1] 


Lf 


^4;zrV 

v  T  y 


(8) 


where  r  is  the  distance  in  meters  between  the  transmitter  and  receiver  antennas  and  A  the 
wavelength,  also  in  meters. 

With  frequency  in  megaHertz  and  the  distance  in  kilometers  d  [km] ,  the  free  space 
loss  in  decibels  is  [11] 

4[dB]  =  32.4  +  20 log d [km]  +  20 log  fc  [MHz],  (9) 


4.  The  Ikegami  Model 

In  this  physical  model,  diffraction  is  included  and  calculated  using  a  single  edge 
approximation  at  the  building  nearest  to  the  receiver.  The  wall  reflection  loss  is  given  as  a 
constant  value.  This  is  shown  in  Figure  5. 


9 


The  path  loss  equation  is  given  as 


I[dB]  =  ZF+4, 


(10) 


where  L¥  is  the  free  space  loss  and  Lcx  is  the  excess  loss  due  to  the  reflection  and  diffrac¬ 
tion. 


Based  on  both  reflected  and  diffracted  rays,  the  equation  for  the  excess  loss  is 
given  as  [1 1] 


4x  [dB]  =  1 0 log  fe  + 1 0  log  (sin  </>)  +  20  log  ( h  -  hm ) 


-20  log  w- 10  log 


'  3  ^ 

1 H — Y 
V  Kj 


-5.8 


(ID 


where  (f)  is  the  angle  between  the  street  and  the  direct  line  of  sight  from  transmitter  to  re¬ 
ceiver.  The  receiver  is  assumed  to  be  in  the  middle  of  the  street  as  shown  in  Figure  6.  In 
Equation  (11),  h  is  the  height  of  the  building  which  causes  the  diffraction,  w  is  the  width 
of  the  street  (between  two  buildings)  where  the  receiver  is  located,  and  L  =  0.25  is  the 
fixed  reflection  loss.  This  model  does  not  consider  the  height  of  the  transmitter  and  the 
building  geometry. 


Figure  6.  Physical  interpretation  of  street  orientation  angle  (j)  (From  Ref.  [8].) 

5.  Walfisch-Bertoni  Model 

In  this  physical  model  [8],  Walfisch  and  Bertoni  accounted  for  the  transmitter 
height  and  building  width,  which  are  not  considered  in  the  Ikegami  model. 
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From  Figure  7,  besides  the  single  diffracted  and  reflected  paths  (labeled  1  and  2), 
the  total  signal  at  the  receiver  may  have  components  from  multiple  rooftop  diffractions 
and  reflections  (labeled  4)  and  building  penetration  (labeled  3).  However,  these  compo¬ 
nents  are  assumed  negligible  [8].  The  path  loss  equation  is  given  as  [8] 

I[dB]  =  32.4  +  20  log  fc  [MHz]  +  20  log  d  [km]  +  Zex  (12) 


where  Lex  is  the  excess  path  loss  due  to  the  diffractions  and  reflections 


Lex  =57.1  +  ^4  +  log/c  +181ogd-181og(/zb -h) 

d 2 


- 1 8  log 


1- 


17(4,-*). 


(13) 


and  A  accounts  for  building  geometry 


A  =  5  log 


f-T 

,2, 


+(*-*, „)! 


9  log  h-l- 20  log 


2{h-hl 

~b~ 


(14) 


Figure  7.  Physical  interpretation  of  Walfisch-Bertoni  model  (From  Ref.  [8].) 


6.  COST  23 1-W alfisch-Ikegami  Model 

This  model  gives  a  better  estimation  of  signal  loss  in  an  urban  environment  be¬ 
cause  it  combines  both  the  Ikegami  and  Wallisch-Bctoni  models.  The  total  loss  equation 
is  given  as  [1 1] 

L[dB]  =  ZI+im-+Zu  (15) 
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where  L¥  is  the  free  space  path  loss,  Lrnsd  is  the  multiple-screen  diffraction  loss  due  to 
multiple  knife-edge  diffraction  to  the  top  of  the  last  building,  and  Lsi  is  for  the  single  dif¬ 
fraction  and  scattering  process  down  to  the  street  level 


Ad  = -16.9  + 101og/c  +10  log 


r(h-hmp 


w 


+  L(<f>) 


where  L  (0)  is  to  account  for  the  street  orientation 


m=\ 


-10  +  0.354 0 
2.5  +  0.075(0-35°) 
4.0-  0.1 14(0-55°) 


for  0°  <  0  <  35° 
for  35°<0<55° 
for  55°  <  0  <  90°. 


(16) 


(17) 


Previously,  in  the  Walfisch-Bcrtoni  model,  the  transmitter  is  above  the  rooftop 
height.  This  model  extends  that  model  to  account  for  a  transmitter  height  below  the  roof¬ 
top  height.  The  loss  equation  given  as  [1 1] 


Ansd  =  4sh  +  K  +  K  lo§  R  +  K  lo§  /c  -  9  l°g  w 

where 

J-181og[l  +  (Ab  -A)]  for  \>h 
[o  for  hb<h, 


(18) 


(19) 


and 


54 

<  54 - 0.8 (\  -h) 

( k—h)R 

54-0.8— - 

{  0.5 

(18 


for  \>h 

for  R  >  0.5  km  and  l\<h 
for  R  <  0.5  km  and  hb  <  h, 

for  hb  >  h 


K  =< 


18-15 


{K-h) 


for  hb  <  h, 


(20) 


(21) 


f 


kf  =-4  +  0.7 


A. 

925 


\ 

1  for  medium-sized  city. 

) 


(22) 
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B.  SHADOW  FADING 

In  an  urban  environment,  there  are  signal  variations  due  to  changes  in  building 
heights,  widths,  shapes,  the  street  width,  and  the  presence  of  street  lamps  and  trees.  The 
distribution  of  such  fading  is  log-normal  because  the  signal,  when  measured  in  decibels, 
has  a  normal  distribution.  The  standard  variation  of  the  log-normal  distribution  for  the 
shadow  fading  is  also  known  as  the  location  variability,  crL  [1 1].  In  mobile  communica¬ 
tions  the  shadow  fading  variation  has  found  to  be  in  the  range  of  6  <  aL  <  12  dB  [12]. 


C.  SMALL-SCALE  FADING 

When  a  transmitter  or  receiver  is  surrounded  by  buildings  or  structures,  they  re¬ 
flect,  scatter  and  diffract  the  transmitted  radiowaves.  The  result  is  that  several  waves  ar¬ 
rive  at  the  receiver  via  different  directions  and  paths.  These  different  radiowaves  can  re¬ 
sult  in  either  signal  strength  addition  or  subtraction  at  the  receiver,  because  the  signals 
have  different  phases,  which  are  independent  and  uniformly  distributed  from  0  to  2 n. 

The  Rayleigh  distribution  is  an  excellent  representation  for  the  signal  strength 
variations  in  NLOS  situations  due  to  small-scale  fading.  It  is  also  known  as  Rayleigh  fad¬ 
ing  [11]. 

In  LOS  situations,  the  received  signal  consists  of  the  Rayleigh  distributed  signal 
plus  a  LOS  signal.  The  resultant  signal  strength  variations  follow  a  Rician  distribution. 
This  Rician  distribution  applies  whenever  one  signal  path  is  much  stronger  than  the  other 
multipath  components  [11]. 


D.  RAY  TRACING  METHOD 

The  previous  sections  have  presented  the  basic  propagation  mechanisms  and  phe¬ 
nomena  such  as  large-scale  path  loss,  the  effect  of  shadowing,  and  finally  the  statistical 
characteristics  of  small-scale  fading.  These  mechanisms  contribute  to  the  path  loss  and 
power  fading  of  radio  signals  in  an  urban  environment. 

The  ability  to  predict  the  fading  is  important  when  planning  a  mobile  radio  system 
to  be  deployed  in  a  particular  environment  or  area.  The  predictions  can  be  included  in  the 
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planning  and  operation  of  such  a  radio  system  to  ensure  efficient  and  effective  communi¬ 
cations  coverage  over  the  desired  area  [8].  Ray  tracing,  which  uses  computational  nu¬ 
merical  methods,  is  an  efficient  way  of  predicting  the  fading. 

Ray  tracing  uses  a  detenninistic  approach.  Rays  are  traced  from  the  transmitter 
and  allowed  to  interact  with  the  environment  yielding  diffracted,  reflected,  and  transmit¬ 
ted  (partially  absorbed  in  the  process),  components  at  a  given  receiver  location.  The 
lengths  of  the  various  ray  paths  give  the  amplitudes  and  phases  of  signal  components  ar¬ 
riving  at  the  receiver.  The  strengths  of  reflected  and  transmitted  rays  are  computed  using 
geometrical  optics  (GO),  and  for  diffracted  rays,  the  geometric  theory  of  diffraction 
(GTD)  [13]. 

1,  Geometrical  Optics 

Geometrical  optics  uses  simple  ray  tracing  to  find  an  approximate  field  strength  at 
the  receiver  using  plane  wave  Fresnel  reflection  and  transmission  formulas  [11].  The  total 
GO  field  is  given  as  [1 1] 

Ni  N, 

EGo  =  E„4e-'v"  +  y  RE  ,4<f'*«  +  £  TE,  A/1'1  (23) 

;=i  j=i 

where  EGO  is  a  2-by-l  column  vector  with  elements  corresponding  to  the  horizontal  and 
vertical  polarized  field  components.  The  other  quantities  in  Equation  (23)  are: 

Nr,Nt,  the  total  numbers  of  reflected  and  transmitted  traced  rays,  respectively, 

rn ,  the  distance  along  nth  ray, 

ki ,  the  wavenumber  associated  with  the  medium  in  which  the  / th  ray  propagates, 

Ai ,  the  spreading  factor  for  the  ith  ray, 

E.  . ,  the  2-by-l  incident  field  vector  at  the  corresponding  transmission  or  reflec¬ 
tion  point, 

R ,  a  2-by-2  reflection  coefficient  matrix,  and 

T  ,  a  2-by-2  transmission  coefficient  matrix. 
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The  spreading  factor,  A.  =  1/r. ,  is  used  for  spherical  waves  and  plane  boundaries. 
Figure  8  shows  some  examples  of  reflection  and  transmission. 


Figure  8.  Propagation  mechanisms  in  real  environment  (From  Ref.  [11].) 

2.  Geometric  Theory  of  Diffraction 

Geometrical  optics,  however,  does  not  include  the  effects  of  ray  diffractions 
caused  by  knife-edges  and  comers  of  structures  in  an  urban  environment.  Keller  devised 
the  geometric  theory  of  diffraction  (GTD)  to  account  for  the  diffractions  [11].  Figure  9 
shows  a  ray  being  diffracted  from  a  wedge. 
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Figure  9.  Edge-diffracted  rays  according  to  GTD  (From  Ref.  [11].) 


The  diffracted  field  is  given  as  [11] 


Erf  =  DE  tAde 


—jkr 


(24) 


where 


f 

nd  II 

,  D  = 

"A 

0  ] 

,  E.  = 

~E*  ' 

H 

1 _ 

0 

Eu._ 

(25) 


and 


Ad=< 


1  /  Vr  for  plane  or  cylindrical  wave  incidence 

\\lri/r{r  +  ri)  t'or  spherical  wave  incidence. 


(26) 


Here,  Ell{  and  E.  denote  the  electric  field  parallel  and  perpendicular  components  relative 
to  the  plane  of  incidence,  respectively,  while  Edn  and  Ed±  are  the  electric  fields  parallel 
and  perpendicular  to  the  plane  of  diffraction.  Also,  Ad  is  the  spreading  factor,  r  is  the 
distance  from  the  source  to  the  diffracting  edge,  r  is  the  distance  from  the  diffracting 
edge  to  the  field  point,  k  is  the  wavenumber,  and  D  and  D  are  the  diffraction  coeffi¬ 
cients. 
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An  extension  of  GTD  is  the  uniform  theory  of  diffraction  (UTD)  which  provides 
more  accurate  prediction  as  it  includes  field  points  in  the  Fresnel-zone-like  region  close 
to  the  shadow  boundary.  Since  UTD  accounts  for  all  points  in  space,  it  is  more  suitable 
for  practical  propagation  applications  [11]. 

Finally,  the  total  field  at  an  observation  point,  P,  can  be  obtained  as  the  sum  of  all 
the  GO  and  GTD  rays  arriving  at  P  [9] 

Et(P)  =  Ego(P)  +  E,(P).  (27) 

In  this  thesis,  a  ray  tracing  engine  (Urbana)  applies  GO  and  GTD  to  model,  simulate  and 
predict  the  performance  of  UAV  data  links  in  an  urban  environment. 

In  this  chapter,  the  concepts  behind  the  mechanisms  of  radiowave  propagation  in 
complex  urban  environment  were  presented.  These  empirical  and  physical  models  can  be 
used  to  describe  the  characteristics  of  radiowave  propagation  in  an  urban  environment. 
Ray  tracing  is  a  useful  tool  for  characterizing  the  radiowave  propagation  channel.  A  ray 
tracing  engine  was  used  for  the  modeling  and  simulation  process  in  this  thesis.  However, 
before  the  simulations  are  presented,  a  summary  of  diversity  techniques  is  given  in  the 
next  chapter. 
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III.  DIVERSITY  TECHNIQUES 


Mobility  is  a  fundamental  feature  of  “untethered”  communications  networks  [14]. 
The  mobility  of  combat  units  is  an  essential  attribute  to  combat  effectiveness.  Thus,  there 
is  great  demand  for  mobile  wireless  communication  systems  in  military  applications.  As 
discussed  in  the  previous  chapter,  radio  propagation  in  urban  areas  is  subjected  to  multi- 
path  and  fading  that  can  severely  degrade  the  quality  of  transmission  and  cause  loss  of  a 
communications  link.  The  challenge  is  how  to  ensure  good  and  reliable  radio  communi¬ 
cations  connectivity  without  limiting  operational  flexibility. 

The  countermeasures  for  fading  can  be  found  in  implementing  diversity  concepts. 
The  basic  theory  of  diversity  is  based  on  the  fact  that  there  is  a  low  chance  of  independ¬ 
ent  fading  signal  paths  experiencing  the  same  fading,  which  means  there  is  a  good  chance 
that  one  of  the  signal  paths  is  not  severely  faded.  Independent  fading  paths  can  be 
achieved  by  separating  the  signal  in  time,  frequency,  space,  or  polarization  [14]. 

A.  TIME  AND  FREQUENCY  DIVERSITY 

Time  diversity  is  achieved  by  repeatedly  transmitting  the  signal  at  a  time  spacing 
equal  to  or  larger  than  the  coherence  time  of  the  channel.  Thus  each  received  signal  is 
subjected  to  independent  fading.  For  example,  in  a  CDMA  RAKE  receiver  implementa¬ 
tion,  the  RAKE  receiver  is  able  to  align  the  multiple  received  signals  in  time  so  that  a  bet¬ 
ter  estimation  of  the  original  signal  may  be  obtained  [15]. 

Frequency  diversity  means  that  the  signal  is  transmitted  using  more  than  one  car¬ 
rier  frequency  with  frequency  spacings  equal  to  or  larger  than  the  coherence  bandwidth 
of  the  channel.  Thus  each  received  signal  is  at  different  carrier  frequency,  and  is  therefore 
subjected  to  independent  fading  [15]. 

Besides  time  and  frequency  diversity,  there  are  other  less  commonly  used  tech¬ 
niques  such  as  spatial  diversity,  cross  polarization  diversity  and  angle  diversity.  These 
were  the  focus  of  the  research  in  this  thesis. 
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B.  SPATIAL  DIVERSITY 

Spatial  diversity  can  be  achieved  with  multiple  combinations  of  transmit  or  re¬ 
ceive  antennas  in  a  radio  system.  Basically,  the  antennas  are  spaced  sufficiently  far  apart 
to  ensure  de-correlation  and  independent  fading  of  the  transmit  or  receive  signals  from 
the  antennas  [10]. 

1.  Diversity  on  Receive 

Independent  fading  paths  in  space  are  achieved  using  a  receiving  antenna  array, 
where  each  element  receives  a  separate  path  signal.  The  antenna  elements  are  required  to 
be  spaced  to  ensure  uncorrelated  received  signals. 

Figure  10  shows  the  relationship  between  omni-directional  antenna  spacing  and 
its  covariance  factor.  For  example,  a  0.4/1  spacing  is  enough  to  ensure  uncorrelated  fad¬ 
ing  between  two  omni-directional  antennas.  A  separation  between  two  antennas  of  one- 
half  to  one  wavelength  is  considered  sufficient  to  obtain  two  nearly  uncorrelated  signals 
[10]. 


Figure  10.  Omni-directional  Antenna  Spatial  Covariance  (After  Ref.  [15].) 

The  larger  the  separation  in  wavelengths,  the  lower  is  the  correlation  between  the 
received  signals.  It  is  important  for  antenna  diversity  that  the  received  signal  be  uncorre¬ 
lated  and  independent  to  ensure  that  there  is  a  chance  of  stronger  signal  on  one  path  if  the 
other  path  is  subjected  to  severe  degradation  [15]. 
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Figure  1 1  shows  the  potential  advantage  of  employing  two  antenna  for  receive  di¬ 
versity.  With  the  antennas  spaced  d  >  X/2  apart  to  ensure  uncorrelated  received  signals, 
the  Gaussian  distributed  received  signals  have  cumulative  distribution  functions  (CDFs) 
that  show  a  diversity  gain,  that  is,  an  improvement  in  the  signal-to-noise  ratio  (SNR).  A 
steep  CDF  is  desirable.  Using  this  advantage,  the  received  signals  can  be  combined  using 
various  techniques  that  can  improve  the  receiver  performance.  The  techniques  are:  (1)  se¬ 
lection  combining  (SC),  (2)  maximal-ratio  combining  (MRC),  and  (3)  equal-gain  com¬ 
bining  (EQC)  [12]. 


Figure  11.  Two  antenna  for  receive  diversity  (After  Ref.  [10].) 

In  the  SC  method  for  two  antennas,  as  shown  in  Figure  12,  the  branch  with  the 
maximum  voltage  SNR  is  selected  as  the  contributing  received  signal.  For  example, 
branch  one  SNRl  is  chosen  if  it  is  larger  than  branch  two  SNR2 .  The  weaker  branch  will 
not  be  used.  In  this  method,  the  mean  SNR  of  the  selected  signal  is  [12] 

N  1 

(SNR)  =  r£-  (28) 

„=1  n 

where  T  is  the  mean  branch  SNR.  The  mean  SNR  improves  logarithmically  with  the 
N  number  of  branches. 


Figure  12.  Two  branch  selection  combining  (SC)  with  equal  noise  power  (After  Ref. 

[16].) 
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In  the  MRC  method  shown  in  Figure  13,  each  branch  signal  is  first  weighted  by 
its  received  instantaneous  voltage  SNR,  r\ / s[~N  and  r2 / ypV~,  respectively.  Thus  the  branch 

with  the  higher  voltage  SNR  is  weighted  more  than  the  branch  with  lower  voltage  SNR. 
The  weighted  signals  are  then  co-phased  and  coherently  added.  This  method  yields  the 
highest  instantaneous  SNR  possible  with  any  linear  combining  technique.  However,  im¬ 
plementation  of  the  MRC  is  expensive  since  the  weights  need  both  amplitude  and  phase 
tracking  of  the  channel  response  [12]. 

SNR]  = 

SNR  = 

2 

Figure  13.  Two  branch  maximal-ratio  combining  (MRC)  with  equal  noise  power  (Af¬ 
ter  Ref.  [16].) 

The  mean  SNR  at  the  output  of  the  combiner  is  [12] 

N 

(SNR)  =  Y,r  =  NT  (29) 

n= 1 

where  T  is  the  mean  branch  SNR.  The  mean  SNR  improves  with  the  number  of 
branches,  N. 

In  the  EQC  method,  the  weights  have  the  same  constant  magnitude,  G,  unlike  the 
MRC,  where  the  weights  are  based  on  the  instantaneous  voltages.  As  shown  in  Figure  14, 
the  weighted  signals  are  then  co-phased  and  coherently  added.  The  noise  power  for  the 
output  voltage  SNR  is  doubled. 


Figure  14. 


Two  branch  equal  gain  combining  (EGC)  with  equal  noise  power  (After 

Ref.  [16].) 
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[12] 


The  closed-form  expression  for  the  mean  SNR  at  the  output  of  the  combiner  is 


(SNR)  = 


i+N-Of 


(30) 


where  T  is  the  mean  branch  SNR.  The  mean  SNR  improvement  is  less  than  the  MRC 
improvement. 

From  Figure  15,  the  plots  of  the  improvement  to  the  mean  SNR  show  that  the 
MRC  and  EQC  have  better  improvement  gain  than  the  SC  method.  For  larger  N,  the 
MRC  is  about  1  dB  better  than  the  EQC  method. 


Figure  15.  Improvement  in  mean  SNR  for  various  combining  techniques  (After  Ref. 

[12].) 


2.  Diversity  on  Transmit 

For  transmit  diversity,  as  shown  in  Figure  16,  the  transmitter  power  is  split  be¬ 
tween  the  antennas.  With  the  antennas  being  spaced  d  >  T/2  apart  to  ensure  uncorrelated 
received  signals,  the  received  signal  strength  (SNR)  which  is  averaged  over  all  angles 
does  not  improve,  but  the  variance  of  the  signal  is  reduced  [10]. 
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Received  Signal  (CDF) 


Figure  16.  Two  antenna  for  transmit  diversity  (After  Ref.  [10].) 


The  advantage  in  SNR  improvement  of  diversity  on  transmit  might  not  be  as  sig¬ 
nificant  as  diversity  on  receive  but  it  has  its  important  use  in  the  multiple  input,  multiple 
output  (MIMO)  concept.  With  multi-element  antenna  (MEA)  arrays  employed  at  both  the 
transmitter  and  receiver  ends,  multipath  in  a  rich  scattering  environment  can  be  employed 
via  space-time  processing  to  achieve  very  high  spectral  efficiencies  [12].  Figure  17  shows 
the  capacity  efficiency  for  MIMO  arrays. 


Achievable  Data  Rates  with  Multiple  Antennas 


Figure  17.  Capacity  efficiency  with  multiple  antennas  (After  Ref.  [17].) 


It  can  be  observed  that  for  a  SNR  per  receiver  antenna  of  about  17  dB  (boxed  lo¬ 
cation)  in  a  30  kHz  channel,  a  single  antenna  transmit-receive  pair  (red  plot)  can  achieve 
a  data  rate  of  about  70  kbps,  while  a  two  antennas  transmit-receive  pair  (brown  plot)  can 
achieve  a  data  rate  of  about  175  kbps.  There  is  an  improvement  of  2.5  times  in  data  rate 
capacity. 
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C.  POLARIZATION  DIVERSITY 

In  addition  to  spatial  diversity,  the  use  of  orthogonal  polarizations  to  exploit  po¬ 
larization  diversity  is  another  approach.  This  approach  provides  for  only  two  diversity 
branches.  Orthogonally  polarized  antennas  ensure  de-correlation  of  the  two  channels 
[15]. 

Figure  18  shows  the  concept  of  polarization  diversity.  A  horizontal  antenna  re¬ 
ceives  only  horizontally  polarized  waves.  Likewise,  a  vertical  antenna  receives  only  ver¬ 
tically  polarized  waves.  In  this  approach,  two  antenna  dipoles  can  be  co-located,  whereby 
one  dipole  is  orientated  vertically  and  the  other  horizontally.  These  orthogonal  antennas 
are  co-located  to  save  space.  This  technique  is  also  known  as  dual-polarization. 


Figure  18.  Polarization  of  waves  generated  by  orthogonal  linearly-polarized  antennas 


D.  ANGLE  DIVERSITY 

Another  diversity  technique  is  the  use  of  directional  antennas  to  achieve  angular 
diversity.  This  approach  is  the  basis  of  space-division  multiple  access  (SDMA).  Angle 
diversity  or  SDMA  can  result  in  an  increase  of  spectral  efficiency  and  capacity  [18]. 

Directional  antennas  have  larger  gains  and  narrower  beamwidths  when  compared 
to  an  omni-directional  antenna.  They  enable  the  transmitter  to  steer  the  antenna  beam  to¬ 
wards  the  location  of  the  intended  receiver.  At  the  same  time,  the  directional  antenna  also 
reduces  reception  by  another  receiver  or  even  interference  by  another  transmitter,  thus 
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providing  a  low  probability  of  detection  (LPD)  system.  With  directional  antennas  having 
sector  patterns,  the  same  frequency  and/or  time  slot  can  be  used  in  each  sector,  as  a  result 
more  users  can  operate  in  the  same  spectrum,  improving  the  system  capacity  [18]. 

Figure  19  shows  the  concept  of  SDMA  with  directional  antennas.  With  the  trans¬ 
mitter  steering  the  antenna  beam  (blue  pattern)  towards  User  1,  it  reduces  reception  and 
interference  by  User  2,  as  it  is  located  in  the  null  of  the  antenna  beam.  Similarly,  when 
the  antenna  beam  (yellow  pattern)  is  steered  towards  User  2,  User  1  has  no  reception  and 
provides  no  interference.  Basically,  SDMA  improves  system  capacity  and  signal  strength 
for  both  transmitter  and  receiver,  at  the  same  time  it  reduces  the  effect  of  interference. 


Figure  19.  Directional  antennas — SDMA  (From  Ref.  [10].) 

In  this  chapter,  the  concepts  behind  the  three  diversity  techniques  of  (1)  spatial 
diversity,  (2)  polarization  diversity,  and  (3)  angle  diversity  or  SDMA,  were  discussed. 
These  diversity  techniques  were  applied  in  the  modeling  and  simulation  of  radiowave 
propagation  in  an  urban  environment.  The  purpose  was  to  investigate  the  effects  of  diver¬ 
sity  on  the  performance  of  urban  radiowave  links.  The  next  chapter  describes  the  soft¬ 
ware  tool  used  in  the  modeling  and  simulation. 
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IV.  MODELING  USING  URBANA 


A.  RAY  TRACING  WITH  URBANA 

The  ray  tracing  software,  Urbana  Wireless  Toolkit,  from  Science  Applications  In¬ 
ternational  Corporation  (SAIC),  was  used  to  model  and  simulate  the  effects  of  diversity 
on  urban  radiowave  links.  It  is  a  computational  electromagnetic  (CEM)  tool  utilizing  a 
ray  tracing  engine  that  combines  physical  optics,  geometric  optics,  and  diffraction  phys¬ 
ics  to  produce  a  3-D  simulation  for  radiowave  propagation  [19]. 

The  deterministic  ray  tracing  techniques  of  Urbana  can  be  used  to  characterize  ra¬ 
diowave  propagation  in  an  urban  environment.  Urbana  uses  ray  tracing  applied  to  CAD 
models  for  buildings,  and  other  major  features.  The  ray  tracing  engine  would  interrogate 
the  3-D  geometry  of  the  environment  to  provide  inputs  to  the  physical  models  for  reflec¬ 
tion,  diffraction  around  buildings  and  terrain,  etc.  [19]. 

For  urban  radiowave  propagation  modeling,  the  input  parameters  for  Urbana  are 
CAD  facet  models  for  the  buildings  and  roads,  including  its  surface  material  (e.g.,  con¬ 
crete,  earth,  glass  and  dielectric  properties).  The  other  inputs  are  the  placement,  strength, 
and  antenna  patterns  of  transmitters  and  receivers.  The  outputs  are  composite  signal 
strength,  angle  of  arrival  and  time  delay  of  each  arriving  signal  ray  at  each  receiver  or  ob¬ 
servation  point  [19]. 


B.  URBAN  CITY  MODEL 

In  this  thesis,  a  small  urban  city  model  adapted  from  an  existing  Urbana  Wireless 
Toolkit  model  was  used.  The  Urbana  model  was  a  replica  of  downtown  Austin,  Texas.  A 
view  of  the  city  model  is  shown  in  Figure  20. 
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Figure  20.  Small  city  model  without  ground  plane 


The  city  model  occupies  a  box  area  of  600  m  by  400  m  by  150  m.  The  tallest 
building  is  at  a  height  of  150  m.  There  are  23  separate  buildings  in  the  model. 


1.  Building  Materials 

Previous  calculations  by  Lock  [20]  and  Pala  [21]  investigated  the  effects  of  vari¬ 
ous  building  materials  on  UAV  data  links.  The  outdoor  signal  levels  are  not  all  that  sensi¬ 
tive  to  the  building  materials  as  long  as  all  of  the  buildings  are  not  completely  perfect 
electric  conducting  (PEC).  Therefore,  since  the  evaluation  of  diversity  techniques  was  the 
goal  of  this  research,  all  the  buildings  were  modeled  with  the  same  material. 

In  Urbana  modeling,  the  23  blocks  of  building  are  made  up  of  3 14  facets  with  ma¬ 
terial  parameters  stated  in  Table  2.  Each  concrete  facet  is  a  penetrable  layered  slab  with 
air  backing. 


COATING 

MATERIAL 

RELATIVE 

PERMITTIVITY 

RE 

PERM 

LATIVE 

IE  ABILITY 

RESISTIVITY 

THICKNESS 

s' 

s" 

m' 

m" 

p  (ohm) 

t(m) 

Concrete 

10.1 

0.5 

1.0 

0 

10JU 

0.3 

Table  2.  City  model  building  material 
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2.  Building  Edges 

Besides  modeling  the  building  surfaces,  the  edges  of  each  building  are  also  in¬ 
cluded  in  Urbana,  as  shown  in  Figure  21.  The  purpose  of  the  edges  is  to  account  for  the 
various  radiowave  diffractions,  which  occur  from  building  and  rooftop  edges. 


Figure  2 1 .  Building  edges  in  the  small  city  model 


3.  Ground  Plane 

A  semi-infinite  ground  plane  is  included  with  the  city  model.  The  city  model  ba¬ 
sically  resides  on  a  semi-infinite  ground  plane  as  shown  in  Figure  22,  which  accounts  for 
ground  reflections.  The  ground  plane  has  s'  =  3  and /u '  =  1 . 


C.  URBANA  USER  INPUT  PARAMETERS 

In  order  for  Urbana  to  model  and  simulate  the  characteristics  of  radiowave  propa¬ 
gation  in  any  environment,  there  are  a  number  of  parameters  that  the  user  needs  to  pro¬ 
vide. 

Table  3  shows  a  list  of  some  user  input  parameters  for  Urbana  and  their  sample 
values.  An  entire  user  input  file  is  listed  in  Appendix  A.  The  settings  in  Table  3  were 
used  throughout  the  simulations.  The  transmitter  and  receiver  antennas  can  be  modeled  as 
dipoles  or  by  using  an  antenna  pattern  data  file  describing  a  typical  radiation  pattern.  The 
modeling  of  the  antennas  is  dependent  on  the  type  of  scenario  to  be  simulated. 


PARAMETERS 

SETTINGS 

Scatterer  model 

city  wgrnd, facet 

Length  unit 

meters 

Frequency 

900  MHz 

Tx  antenna 

dipole  or  antenna  pattern  data 

Rx  antenna 

dipole  or  antenna  pattern  data 

Ray  tracing  computation 

geometrical  optics  (GO) 

Edge  diffraction 

uniform  theory  of  diffraction  (UTD) 

Edge  model 

city,  edge 

Burst  ray  angular  spacing 

2° 

Maximum  ray  bounces 

7 

Coating  material 

from  Table  2 

Advance  features 

off 

Table  3.  Urbana  user  input  parameters 


D.  MODELING  OF  ANTENNAS 

The  characteristics  of  the  transmitting  and  receiving  antennas  are  essential  pa¬ 
rameters  for  the  Urbana  processing.  Basically,  different  types  of  antennas  affect  radio¬ 
wave  propagation  differently  through  their  radiation  patterns.  There  were  two  basic  types 
of  antennas  modeled  in  this  thesis,  discussed  in  the  following  sections. 
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1.  Dipole  Antenna 

TheT/2  dipole  is  the  basic  antenna  element  used  throughout  this  thesis.  The  radia¬ 
tion  pattern  of  aT/2  z-directed  dipole  at  an  operating  frequency  of  900  MHz  in  the 
xz  plane  is  shown  in  Figure  23.  It  has  an  omni-directional  pattern  in  xy  plane. 

Relative  Voltage  Plot 
90  1 


Figure  23.  Half-wave  dipole  radiation  pattern  inxz  plane 

2.  Directional  Antenna 

A  directional  antenna  pattern  was  created  for  use  in  Urbana.  Figure  24  shows  a 
plot  of  the  radiation  pattern.  It  has  a  3-dB  half-power  beamwidth  (HPBW)  of  10°  at  an 
operating  frequency  of  900  MHz.  The  radiation  into  the  rear  hemisphere  was  set  to  zero, 
because  it  was  assumed  that  this  pattern  is  due  to  an  array  of  elements  with  a  cavity  or 
ground  plane  backing.  The  directivity  of  the  directional  antenna  is  about  25  dB. 
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Figure  24.  Directional  antenna  radiation  pattern  with  FIPBW=  10° 


In  Urbana,  an  antenna  pattern  file  containing  the  linear  field  strength  components 
of  the  above  directional  antenna  was  created.  The  pattern  file  is  generated  using  the  Mat- 
lab  codes  shown  in  Appendix  B.  The  pattern  is  nonnalized  to  the  directivity,  which  is 
automatically  computed  by  Urbana.  Figure  25  shows  the  free  space  radiation  pattern 
simulated  in  Urbana  with  an  excitation  of  1  W,  and  pointed  in  the  direction  of  an¬ 
gles  ^  =  0°  and  6  =  90°,  respectively.  The  location  of  the  antenna  is  at  (x,y,z)  =  (-10  m, 
-55  m,  - 150  m).  The  observation  plane  is  at  a  height  of  z  =  4  m. 
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Figure  25. 


Signal  coverage  in  the  forward  hemisphere  of  the  directional  antenna  at 
location (-10  m,  -55  m,  150  m) 
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E.  MODELING  OF  DIVERSITY  TECHNIQUES  FOR  UAV 

COMMUNICATIONS  LINKS 

With  an  appreciation  of  the  UAV  operations  and  the  critical  need  for  effective 
tactical  communications,  we  now  explore  methods  of  implementing  various  diversity 
techniques  to  improve  the  effectiveness  of  UAV  communication  links  in  MOUT.  This 
also  sets  the  framework  for  simulations  later  on. 

1.  Spatial  Diversity 

In  order  to  investigate  the  benefit  of  antenna  spatial  diversity,  two  antennas  were 
considered  in  the  simulations.  Further  improvement  is  achievable  with  four,  eight  or  even 
more  antennas,  but  at  the  expense  of  increased  antenna  area  (“real  estate”). 

a.  Spatial  Diversity  on  Receive 

In  the  Urbana  modeling  process,  a  vertical  A/2  dipole  was  modeled  as 
transmitting  from  a  UAV  hovering  150  m  above  the  city  center,  as  shown  in  Figure  26. 
The  UAV  structure  itself  was  not  considered.  It  was  assumed  that  the  dipole  type  pattern 
is  achieved  when  the  antenna  is  installed  on  the  UAV. 


Figure  26.  Vertical  T/2  dipole  transmitting  at  location  (-10  m,-55  m,  150  m) 

On  the  ground,  observation  points  were  located  at  a  height  of  1 .5  m  above 
the  ground  plane  and  covered  the  600  m  by  400  m  small  city  model  with  cell  size  of  2  m. 
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Some  of  these  observation  points  were  designated  as  receivers  with  2/2  dipoles  for  de¬ 
tailed  analysis. 

After  the  simulations,  spatial  diversity  on  receive  was  assessed  by  analyz¬ 
ing  the  signal  strength  obtained  at  the  receivers.  A  two-antenna  diversity  approach  on  re¬ 
ceive  can  be  evaluated  by  analyzing  a  pair  of  the/1/2  dipoles  which  are  placed  2  m  apart. 
The  dipole  separation  is  more  than  102,  apart  at  the  operating  frequency  of  900  MHz. 
This  ensures  good  de-correlation  between  the  received  signals.  Additional  simulations 
and  analysis  with  different  orientations  of  the  transmitting  and  receiving  dipoles  were 
also  investigated. 

b.  Spatial  Diversity  on  Transmit 

In  modeling  the  spatial  diversity  on  transmit,  a  pair  of  2/2  dipoles  were 
modeled  on  the  transmitting  UAV.  The  dipoles  were  placed  at 2/2  apart  (0. 167  m  at  900 
MHz)  to  ensure  sufficient  de-correlation  of  the  transmit  signals.  Figure  27  shows  the  spa¬ 
tial  diversity  model  on  transmit.  In  addition,  the  total  excitation  power  for  the  pair  of  di¬ 
poles  was  1  W.  Thus  each  dipole  transmits  0.5  W.  As  a  result,  there  is  a  common  basis 
for  comparison  with  a  single  dipole  transmitter  excited  by  1  W. 


Figure  27.  Spatial  diversity  on  transmit:  two  vertical  dipoles  at  locations 
(-10  m,-55  m,  150  m)  and  (-10.167  m,-55  m,  150  m) 
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On  the  ground,  observation  points  were  located  at  a  height  of  1 .5  m  above 
the  ground  plane  and  covered  the  600  m  by  400  m  small  city  model  with  cell  size  of  2  m. 
Some  of  these  observation  points  are  designated  as  receivers  with  2/2  dipoles  for  detailed 
analysis.  The  advantage  of  having  spatial  diversity  on  both  transmit  and  receive  was  in¬ 
vestigated  by  analyzing  the  total  received  signal  strength. 

2.  Polarization  Diversity 

In  modeling  polarization  diversity,  there  were  two  orthogonal  transmit  dipoles  in 
co-location.  One  antenna  was  vertically  polarized  while  the  other  was  horizontally 
polarized,  as  shown  in  Figure  28.  This  is  also  referred  to  as  dual-polarization  diversity. 


Figure  28.  Dual-polarization  diversity  on  transmit  at  location  (-10  m,-55  m,  150  m) 

On  the  ground,  observation  points  were  located  at  a  height  of  1.5  m  above  the 
ground  plane  and  covered  the  600  m  by  400  m  small  city  model  with  cell  size  of  2  m. 
Some  of  these  observation  points  are  designated  as  receivers  with  2/2  dipoles  for  detailed 
analysis.  The  receiver  antennas  are  A/2  dipoles  orientated  either  vertically,  horizontally, 
or  both.  Analysis  was  conducted  by  comparing  the  performance  of  polarization  diversity 
to  the  spatial  diversity  approach  described  previously. 

3.  Angle  Diversity 

In  order  to  demonstrate  angle  diversity,  a  vertical  A/2  dipole  was  used  as  the 

transmitting  antenna  on  the  UAV.  A  receiver  with  a  directional  antenna  pattern  (shown  in 
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Figure  25)  was  located  on  the  ground.  For  angle  diversity  analysis,  a  few  transmitting  di¬ 
pole  locations  were  selected  and  transmissions  from  the  UAV  simulated.  This  represents 
samples  of  transmissions  along  a  UAV  flight  path.  For  the  receiving  antenna  on  the 
ground,  a  shadow  location  was  chosen  at  a  point  with  deep  fading.  The  receiving  direc¬ 
tional  antenna  was  steered  towards  the  direction  of  the  strongest  incoming  signal,  which 
is  generally  not  in  the  unobstructed  line-of-sight  direction.  In  addition,  a  comparative 
analysis  was  done  for  all  of  the  diversity  techniques  described  previously,  to  assess  their 
relative  perfonnance  in  the  urban  environment. 

In  this  chapter,  the  framework  for  the  simulation  process  was  defined.  The  city 
model  and  the  different  types  of  antennas  to  be  used  were  created  and  modeled.  Also,  the 
scenarios  for  applying  the  diversity  techniques  to  the  UAV  data  links  over  the  urban  city 
were  also  described.  The  next  chapter  presents  the  simulation  results. 
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V.  SIMULATIONS  AND  ANALYSIS 


A.  SIMULATIONS  OF  SPATIAL  DIVERSITY 

Antenna  spatial  diversity  was  simulated  by  using  two  A/2  dipoles.  Antenna  spatial 
diversity  can  be  simulated  on  either  receive,  or  transmit,  or  both.  On  transmit,  antennas 
were  placed  A/2  apart  (0. 167  m  at  900  MHz)  to  ensure  de-correlation  of  the  receive  sig¬ 
nals. 


1.  Spatial  Diversity  on  Receive 

A  UAV  transmitter  was  modeled  as  a  single  dipole  transmitting  over  the  city 
model.  The  observation  points  for  receiving  signals  were  positioned  throughout  the  city 
area  of  600  m  by  400  m  with  cell  size  of  2  m,  and  at  a  height  of  1 .5  m.  A  single  A/ 2  di¬ 
pole  antenna  was  modeled  as  receiver  in  some  of  the  observation  cells  for  analysis.  Also, 
vertical  (V)  and  horizontal  (H)  orientations  of  the  transmitting  and  receiving  dipoles  were 
considered.  Note  that  vertical  dipoles  are  parallel  to  the  z  axis,  while  horizontal  dipoles 
are  parallel  to  they  axis.  Table  4  shows  the  parameters  used  in  the  simulation. 


PARAMETERS 

SETTINGS 

Operating  Frequency 

900  MHz 

Type  of  Tx  antenna 

dipole 

Length  of  Tx  antenna 

0.1667  m 

Location  of  Tx  antenna 

(x,y,z)  =  (-10  m,-55  m,  150  m) 

Tx  power/phasing 

1Z0°  w 

Orientations  of  Tx  antenna 

parallel  to  the  z  axis  /or 
parallel  to  the  y  axis 

Observation  points  (used  for  signal 
coverage  plot  computation) 

600  m  by  400  m  with  cell  size  of  2  m, 
height  of  1.5  m 

Type  of  Rx  antenna 

dipole 

Length  of  Rx  antenna 

0.1667  m 

Locations  of  Rx  antenna 

(145  m,200  m,  1.5  m)  to 

(145  m,-200  m,  1.5  m)  with  2  m  resolution 

Orientations  of  Rx  antenna 

parallel  to  thez  axis  /or 
parallel  to  the  y  axis 

Table  4.  Urbana  input  parameters  for  simulation  of  spatial  diversity  on  receive 
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The  total  signal  strength  coverage  plots  are  shown  in  Figures  29  through  3 1 .  Ur- 
bana  computes  the  total  electric  field  strength  of  all  arriving  signal  rays  at  each  observa¬ 
tion  point  to  plot  out  the  total  signal  strength  in  dBm. 

Figures  29  and  30  illustrate  the  signal  coverage  over  the  city  for  a  vertical  trans¬ 
mitting  dipole.  Figure  29  shows  the  3-D  view  with  building  edges.  The  location  of  the 
transmitting  antenna  is  marked  by  the  symbol  “X”  in  Figure  30.  In  the  urban  environ¬ 
ment,  the  presence  of  buildings  results  in  multipath.  Thus,  there  are  locations  where  fad¬ 
ing  is  dominant  due  to  the  various  mechanisms  discussed  in  Chapter  II.  This  can  be  ob¬ 
served  in  the  variation  of  the  signal  strengths  which  are  represented  by  the  different  color 
codes.  The  strongest  received  signal  is  about -21  dBm  and  weakest  is  -82  dBm. 


Figure  29.  Signal  coverage  (3-D  view  with  building  edges)  for  a  vertical  transmit  di¬ 
pole  at  location  (-10  m,-55  m,  150  m) 
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Figure  30.  Signal  coverage  for  a  vertical  transmit  dipole  at  location 

(-10  m,-55  m,  150  m) 


The  total  signal  strength  coverage  with  a  horizontal  transmitting  dipole,  shown  in 
Figure  31,  ranges  from  a  strongest  received  signal  of  about -14  dBm  to  a  weakest  of 
about  -80  dBm.  The  stronger  signals  have  larger  values  than  those  for  the  vertical 
transmitting  dipole.  However,  the  stronger  signals  are  concentric  near  and  under  the 
transmitter  location  because  the  horizontal  dipole  has  a  downward  radiation  at  the  an¬ 
tenna  location.  A  vertical  dipole  has  a  null  pointing  downward.  Horizontal  radiation  is 
not  as  far  reaching  as  vertical  radiation  most  likely  due  to  the  large  number  of  vertical 
edges,  which  enhances  edge  diffraction. 
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y 

Figure  3 1 .  Signal  coverage  for  horizontal  transmit  dipole  at  location 

(-10  m,-55  m,  150  m) 


In  order  to  establish  a  common  basis  for  comparing  the  different  diversity  tech¬ 
niques,  observation  points  from  (145  m,  200  m,  1.5  m)  to  (145  m, -200  m,  1.5  m)  at  2  m 
resolution  on  the  y  axis  were  used.  This  is  called  a  line  path  in  Urbana,  and  is  shown  by 
the  dashed  line  in  Figure  30.  There  are  200  observation  points  on  this  path.  It  was  chosen 
because  the  reception  is  strong  at  some  points  and  experiences  deep  fading  at  others. 

Each  observation  point  was  modeled  with  a  dipole  and  receiver. 

Figure  32  compares  the  difference  between  a  vertical  dipole  transmitter  and  re¬ 
ceiver  pair  (VTVR)  to  a  horizontal  dipole  transmitter  and  receiver  pair  (HTFIR).  It  is  ob¬ 
served  that  at  locations  near  to  the  transmitter  (observation  points  50  to  125),  the  horizon¬ 
tal  dipole  pair  has  higher  received  signal  strength  as  compared  to  the  vertical  dipole  pair. 
At  locations  farther  away  (observation  points  0  to  50  and  125  to  200),  the  vertical  dipole 
pair  has  better  performance.  There  are  some  deep  fades  encountered  in  the  horizontal 
case  due  to  blockage  by  buildings.  The  horizontal  transmitting  dipole  has  strong  radiation 
toward  the  ground  below  the  antenna  location.  Since  the  pattern  nulls  are  pointed  outward 
toward  the  horizon,  the  horizontal  radiation  is  not  as  far  reaching  as  a  vertical  transmit¬ 
ting  dipole  in  these  directions. 
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Observation  points 

Figure  32.  Comparison  between  vertical  and  horizontal  polarized  links  from  the  UAV 

to  the  ground 

Figures  33  and  34  show  the  results  for  two  receiving  dipoles  in  each  case.  The 
signal  power  for  receive  spatial  diversity  was  computed  using  the  EGC  diversity  tech¬ 
nique  (discussed  in  Chapter  III).  For  the  EGC  diversity  technique,  the  received  signal 
from  each  dipole  was  normalized  by  a  constant  gain,  and  then  co-phased  and  coherently 
added  to  give  a  better  signal  strength. 


Figure  33.  Spatial  diversity  on  receive,  vertical  case:  diversity  vs.  no  diversity 
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Figure  34.  Spatial  diversity  on  receive,  horizontal  case:  diversity  vs.  no  diversity 

The  improvement  using  spatial  diversity  is  obvious.  The  improvement  factor  for 
diversity  on  receive  varies  depending  on  location.  At  a  location  of  deep  fading,  the  im¬ 
provement  in  received  signal  strength  for  both  vertical  and  horizontal  spatial  diversity  is 
about  6  dB,  as  compared  to  no  diversity.  It  is  to  be  noted  that  with  EGC  at  the  receivers, 
noise  power  would  also  be  doubled  (3  dB),  so  the  SNR  improvement  would  only  be  3  dB 
In  both  cases,  spatial  diversity  has  also  reduced  the  number  of  drastic  drops  in  signal 
strength  (fading),  which  occur  regularly  along  the  observation  points  in  the  no  diversity 
scenario. 

Figure  35  shows  the  comparison  between  a  vertical  dipole  transmitter  and  two 
vertical  dipoles  receiver  on  diversity  (VTRD),  to  a  horizontal  dipole  transmitter  and  two 
horizontal  dipoles  receiver  on  diversity  (HTRD).  The  major  difference  is  at  observation 
points  near  to  the  transmitting  location  (50  to  125).  The  horizontal  case  has  better  re¬ 
ceived  signal  power.  However,  at  locations  far  away  from  the  transmitter  (125  to  200), 
the  vertical  case  has  some  improvement  over  the  horizontal  dipoles. 
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Figure  35.  Spatial  diversity  on  receive:  comparison  between  vertical  and  horizontal 

cases 

2.  Spatial  Diversity  on  Transmit 

In  this  scenario,  the  transmitter  uses  antenna  spatial  diversity.  At  the  transmitter, 
two  vertical  A/2  dipoles  were  placed  atA/2  apart  (0.167  m  at  900  MHz),  with  each  being 

excited  by  a  zero-phase  half  unit  power.  Each  dipole  transmitted  independently.  Also, 
vertical  and  horizontal  orientations  of  the  transmitting  and  receiving  dipoles  were  simu¬ 
lated.  The  parameters  used  in  the  simulation  are  shown  in  Table  5. 
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PARAMETERS 

SETTINGS 

Operating  Frequency 

900  MHz 

Type  of  Tx  antenna 

2  dipoles 

Length  of  Tx  antenna 

0.1667  m 

Locations  of  Tx  antennas 

(-10  m,-55  m,  150  m)  and 
(-10.167  m,-55  m,  150  m) 

Tx  power/phasing 

0.5Z00  W  and  0.5Z00  W 

Orientations  of  Tx  antennas 

parallel  to  the  z  axis  /or 
parallel  to  the  y  axis 

Observation  points  (used  for  signal 

600  m  by  400  m  with  cell  size  of  2  m, 

coverage  plot  computation) 

height  of  1.5  m 

Type  of  Rx  antenna 

dipole 

Length  of  Rx  antenna 

0.1667  m 

Locations  of  Rx  antennas 

(145  m,200  m,  1.5  m)  to  (145  m,-200  m,  1.5  m) 
with  2  m  resolution 

Orientations  of  Rx  antennas 

parallel  to  the  z  axis  /or 
parallel  to  the  y  axis 

Table  5.  Urbana  input  parameters  for  simulation  of  spatial  diversity  on  transmit 

The  total  signal  strength  coverage  plots  are  shown  in  Figures  36  and  37.  Figure  36 
provides  the  signal  coverage  over  the  city  for  the  vertical  transmitting  dipoles.  Consider¬ 
ing  all  the  observation  points,  the  strongest  received  signal  is -18  dBm  and  the  weakest  is 
-87  dBm .  By  comparing  this  with  the  level  obtained  with  a  single  transmitting  dipole, 
the  strongest  signal  improves  by  about  3  dB.  The  signal  strength  coverage  with  the  hori¬ 
zontal  transmitting  dipoles,  as  shown  in  Figure  37,  ranges  from  the  strongest  received 
signal  of  about- 13  dBm  to  the  weakest  of  about -85  dBm. 
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Figure  36.  Signal  coverage  for  two  vertical  transmit  dipoles  at  locations 
(-10  m,-55  m,  150  m)  and(-10.167  m,-55  m,  150  m) 
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Signal  coverage  for  two  horizontal  transmit  dipoles  at  locations 
(-10  m,-55  m,  150  m)  and(-10.167  m,-55  m,  150  m) 


Figures  38  and  39  show  the  comparison  of  the  received  signal  power  when  using 
spatial  diversity  on  transmit  for  vertical  and  horizontal  cases,  respectively.  The  two 
curves,  shown  for  each  case,  are  for  diversity  on  receive  (red  plot)  and  no  receive  diver- 
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sity  (blue  plot).  The  signal  powers  for  the  receive  diversity  are  computed  assuming  EGC 
technique  at  the  receivers.  The  perfonnances  for  both  cases  are  similar  to  that  of  having 
no  transmitter  diversity  as  discussed  in  the  following  paragraph. 
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Figure  38.  Spatial  diversity  on  transmit,  vertical  case:  comparison  of  diversity  on  re¬ 
ceive  vs.  no  diversity 
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Figure  39.  Spatial  diversity  on  transmit,  horizontal  case:  comparison  of  diversity  on 

receive  vs.  no  diversity 
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Figures  40  and  41  show  a  comparison  of  diversity  on  transmit  (red  plot)  and  no 
transmit  diversity  (blue  plot)  for  vertical  and  horizontal  cases,  respectively.  It  is  observed 
that  transmit  diversity  has  only  slight  improvement  over  the  no  diversity  case.  As  ex¬ 
plained  in  Chapter  III  previously,  spatial  diversity  on  transmit  has  its  primary  advantage 
in  MIMO  where  high  spectral  efficiency  can  be  achieved. 


Figure  40.  Spatial  diversity  on  receive,  vertical  case:  comparison  between  diversity 

on  transmit  and  no  diversity 


Figure  41.  Spatial  diversity  on  receive,  horizontal  case:  comparison  between  diversity 

on  transmit  and  no  diversity 
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B.  SIMULATIONS  OF  POLARIZATION  DIVERSITY 

Antenna  dual-polarization  diversity  was  simulated  with  two  orthogonally  polar¬ 
ized  co-located  A/2  dipoles,  with  one  orientated  horizontally,  the  other  vertically.  This 
method  is  advantageous  when  area  is  limited  on  the  platfonn  and  spatially  dispersed  an¬ 
tennas  cannot  be  installed.  A  few  scenarios  using  various  antenna  configurations  on  both 
transmit  and  receive  were  investigated. 


1.  Polarization  Diversity  on  Transmit  with  Horizontal  Dipole  on  Receive 

In  this  scenario,  the  UAV  transmitter  antenna  was  modeled  as  two  orthogonally 
polarized  A/2  dipoles.  Each  dipole  was  excited  by  a  zero-phase  half  unit  power.  The  re¬ 
ceiver  antenna  is  a  single  horizontally  polarized  A/2  dipole.  Table  6  shows  the  parameters 
used  in  the  simulation. 


PARAMETERS 

SETTINGS 

Operating  Frequency 

900  MHz 

Type  of  Tx  antenna 

2  dipoles 

Fength  of  Tx  antenna 

0.1667  m 

Focations  of  Tx  antenna 

(-10  m,-55  m,  150  m)  and 
(-10  m,-55  m,  150  m) 

Tx  power/phasing 

0.5Z00  W  and  0.5Z00  W 

Orientations  of  Tx  antenna 

parallel  to  the  z  axis  and 
parallel  to  the  y  axis 

Observation  points  (used  for  signal 
coverage  plot  computation) 

600  m  by  400  m  with  cell  size  of  2  m, 
height  of  1 .5  m 

Type  of  Rx  antenna 

dipole 

Fength  of  Rx  antenna 

0.1667  m 

Locations  of  Rx  antenna 

(145  m,200  m,  1.5  m)  to  (145  m,-200  m,  1.5  m) 
with  2  m  resolution 

Orientation  of  Rx  antenna 

parallel  to  the  y  axis 

Table  6.  Urbana  input  parameters  for  simulation  of  dual-polarization  diversity  on 

transmit  with  horizontal  dipole  on  receive 

The  total  signal  strength  plot  is  shown  in  Figure  42.  It  is  observed  that  strong  sig¬ 
nals  are  concentric  around  the  transmitter  location.  The  strongest  received  signal  strength 
is  about -16  dBm  while  the  weakest  is -81  dBm. 
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(300,200) 


16  dBm 


-81  dBm 


Figure  42.  Signal  coverage  plot  for  dual-polarized  transmit  dipoles  at  location 

(-10  m,-55  m,  150  m) 

For  further  analysis,  a  horizontal  dipole  on  receive  was  simulated  over  the  line 
path  of  200  observation  points.  Figure  43  shows  the  received  signal  strength  on  the  line 
path  from  (145  m,  200  m,  1.5  m)  to  (145  m, -200  m,  1.5  m)  at  2  m  increments. 
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Figure  43.  Dual-polarized  transmit  antenna  with  a  horizontal  dipole  on  receive 
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The  received  total  power  is  contributed  by  both  the  vertical  element  (VT)  and 
horizontal  element  (HT)  of  the  dual-polarized  transmitter  antenna.  It  is  observed  that  the 
main  contribution  to  the  total  signal  strength  is  from  the  horizontally  polarized  element  of 
the  dual-polarized  transmitter  antenna.  This  is  expected  as  the  receivers  are  horizontally 
polarized  antennas. 


2.  Polarization  Diversity  on  Transmit  with  Vertical  Dipole  on  Receive 

In  this  scenario,  the  receiver  antennas  were  configured  to  be  vertically  polarized 
X/2  dipoles.  Table  7  shows  the  parameters  used  in  the  simulation. 


PARAMETERS 

SETTINGS 

Operating  Frequency 

900  MHz 

Type  of  Tx  antenna 

2  dipoles 

Length  of  Tx  antenna 

0.1667  m 

Locations  of  Tx  antenna 

(-10  m,-55  m,  150  m)  and  (-10  m,-55  m,  150  m) 

Tx  power/phasing 

0.5Z00  W  and  0.5Z00  W 

Orientations  of  Tx  antenna 

parallel  to  the  z  axis  and 
parallel  to  the  y  axis 

Type  of  Rx  antenna 

dipole 

Length  of  Rx  antenna 

0.1667  m 

Locations  of  Rx  antenna 

(145  m,  200  m,  1.5  m)  to  (145  m,-200  m,  1.5  m) 
with  2  m  resolution 

Orientation  of  Rx  antenna 

parallel  to  the  y  axis 

Table  7.  Urbana  input  parameters  for  simulation  of  dual-polarization  diversity  on 

transmit  and  a  vertical  dipole  on  receive 

Vertical  receiving  dipoles  were  simulated  over  the  line  path  from  (145  m,  200  m, 
1.5  m)  to  (145  m,  -200  m,  1.5  m)  at  2  m  increments.  The  results  are  shown  in  Figure  44. 
The  received  total  power  is  contributed  by  both  the  vertical  (VT)  and  horizontal  (HT) 
elements  of  the  dual-polarized  transmitter  antenna.  The  main  contribution  to  the  total  sig¬ 
nal  strength  is  from  the  vertically  polarized  element  (VT)  of  the  dual-polarized  transmit¬ 
ter  antenna.  This  is  expected  because  the  receivers  are  vertically  polarized  dipoles. 
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Figure  44.  Dual-polarized  transmit  antenna  with  a  vertical  dipole  on  receive 


3.  Polarization  Diversity  on  both  Transmit  and  Receive 

In  this  scenario,  dual-polarization  diversity  was  simulated  for  both  the  transmitter 
and  receiver  antennas.  The  dual-polarized  receiver  dipoles  were  simulated  over  the  line 
path.  Table  8  shows  the  parameters  used. 


PARAMETERS 

SETTINGS 

Operating  Frequency 

900  MHz 

Type  of  Tx  antenna 

2  dipoles 

Length  of  Tx  antenna 

0.1667  m 

Locations  of  Tx  antenna 

(-10  m,-55  m,  150  m)  and 
(-10  m,-55  m,  150  m) 

Tx  power/phasing 

0.5Z00  W  and  0.5Z00  W 

Orientations  of  Tx  antenna 

parallel  to  the  z  axis  and 
parallel  to  the  y  axis 

Type  of  Rx  antenna 

2  dipoles 

Length  of  Rx  antenna 

0.1667  m 

Locations  of  Rx  antenna 

(145  m,  200  m,  1.5  m)  to  (145  m,-200  m,  1.5  m) 
with  2  m  resolution 

Orientations  of  Rx  antenna 

parallel  to  the  z  axis  and 
parallel  to  the  y  axis 

Table  8.  Urbana  input  parameters  for  simulation  of  dual-polarization  on  transmit 

and  receive 
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The  results  of  this  simulation  were  combined  with  the  results  of  the  previous  sce¬ 
narios.  A  comparative  analysis  was  done  to  assess  the  performance  of  various  combina¬ 
tions  of  polarization  diversity.  These  are  discussed  in  the  following  paragraphs. 

Figure  45  compares  between  dual-polarized,  vertical,  and  horizontal  dipoles  on 
receive,  with  the  transmitter  on  dual-polarized  dipoles.  EGC  technique  was  assumed  to  be 
implemented  at  the  receivers.  It  is  observed  that  a  horizontal  dipole  on  receive  (FIT)  has 
the  advantage  over  a  vertical  dipole  (VT)  at  observation  points  near  to  the  transmitter  lo¬ 
cation  (50  to  125).  However,  for  dual-polarized  dipoles  on  receive  (DPR),  the  received 
signal  power  is  improved  relative  to  a  single  polarization.  For  example,  at  the  deep  fading 
location,  the  dual-polarized  receiver  has  about  a  5-dB  signal  strength  advantage  over  the 
horizontally  polarized  receiver.  However,  it  is  to  be  noted  that  the  SNR  improvement 
would  be  2  dB  as  noise  power  is  doubled  (3  dB)  in  EGC. 
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Figure  45.  Dual-polarization  diversity  on  transmit:  comparison  between  dual- 
polarized,  vertical,  and  horizontal  dipoles,  on  receive 

In  Figure  46,  the  case  of  dual-polarized  dipoles  on  transmit  and  receive  (DPTR)  is 
compared  to  a  vertical  dipole  on  transmit  and  receive  (VTR),  and  a  vertical  dipole  on 
transmit  and  dual-polarized  dipoles  on  receive  (VTDPR).  The  plot  shows  that  the  dual- 
polarized  dipoles  on  transmit  and  receive  (red  plot)  has  better  received  signal  power  at 
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most  locations.  However,  at  the  deep  fading  location,  there  is  no  advantage  of  the  dual- 
polarized  dipoles  on  transmit  over  the  other  two  configurations.  In  addition,  at  observa¬ 
tion  points  (125  to  200),  it  has  a  number  of  drastic  drops  of  signal  strength  (fading). 
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Figure  46.  Comparison  between  dual-polarized  and  vertical  dipole  on  transmit  and 

receive 

Figure  47  compares  the  case  where  dual-polarized  dipoles  on  transmit  and  receive 
(DPTR),  against  horizontally  (DPTHDR)  and  vertically  (DPTVDR)  polarized  spatial  di¬ 
versity  on  receive.  At  the  deep  fading  location,  the  vertical  case  has  a  slight  improved 
signal  power  of  about  3  dB  as  compared  to  dual-polarized  diversity  on  receive.  The  hori¬ 
zontally  polarized  spatial  diversity  gives  improved  signal  strength  at  observation  points 
near  to  the  transmitter  (50  to  125).  For  dual-polarization  on  transmit,  it  is  obvious  that 
spatial  diversity  on  receive  has  an  advantage  over  dual-polarization  on  receive. 
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Figure  47.  Dual-polarization  diversity  on  transmit:  comparison  between  dual- 
polarized  dipoles  and  spatial  diversity  on  receive 

Finally,  a  comparison  was  done  for  spatial  diversity  and  polarization  diversity  on 
both  transmit  and  receive.  In  Figure  48,  the  spatial  diversity  (red  plot)  performs  better 
than  the  dual-polarization  diversity  (blue  plot).  For  example,  at  the  deep  fading  point,  the 
spatial  diversity  technique  gives  about  7  dB  improvement.  At  observation  points  near  to 
the  transmitter  (50  to  125),  both  give  comparable  signal  strengths.  In  addition,  the  num¬ 
bers  of  drastic  signal  drops  (fading)  have  reduced  significantly  for  the  spatial  diversity  as 
compared  to  the  polarization  diversity  case. 
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Figure  48.  Comparison  between  dual-polarization  diversity  and  spatial  diversity  on 

both  transmit  and  receive 

C.  SIMULATIONS  OF  ANGLE  DIVERSITY 

In  the  simulation  of  angle  diversity,  a  directional  antenna  was  modeled  at  the  re¬ 
ceiver  while  a/t/2  dipole  was  modeled  on  the  transmitting  UAV.  The  characteristic  of  the 
directional  antenna  was  discussed  in  Chapter  IV.  The  radiation  patterns  are  shown  in  Fig¬ 
ures  24  and  25.  In  order  to  investigate  the  effects  of  angle  diversity,  the  transmitter  loca¬ 
tions  were  taken  at  discrete  points  along  a  flight  path. 

1.  Angle  Diversity  on  Receive 

Table  9  shows  the  parameters  used  for  simulation  of  a  single  dipole  on  transmit 
with  angle  diversity  on  receive.  There  are  five  transmitter  locations.  The  location  of  the 
receiver  is  chosen  to  be  at  (145  m,  52  m,  1.5  m).  This  location  was  chosen  from  the  line 
path  of  observation  points  that  were  used  in  the  previous  diversity  analyses.  The  purpose 
was  to  assess  the  performance  of  angle  diversity  over  the  other  techniques,  which  did 
very  little  to  mitigate  the  deep  fading. 
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PARAMETERS 

SETTINGS 

Operating  Frequency 

900  MHz 

Type  of  Tx  antenna 

dipole 

Length  of  Tx  antenna 

0.1667  m 

Locations  of  Tx  antenna 

Tx  #1  at  (-80  m,  45  m,  150  m) 

Tx  #2  at  (-45  m, -5  m,  150  m) 

Tx  #3  at  (-10  m, -55  m,  150  m) 

Tx  #4  at  (25  m, -105  m,  150  m) 

Tx  #5  at  (60  m,-155  m,  150  m) 

Tx  power/phasing 

1Z0°  W 

Orientation  of  Tx  antenna 

parallel  to  the  z  axis 

Observation  points  (used  for  signal 
coverage  plot  computation) 

600  m  by  400  m  by  1 .5  m 
with  2  m  resolution 

Type  of  Rx  antenna 

directional  (See  Figure  24) 

Location  of  Rx  antenna 

(145  m,  52  m,  1.5  m) 

Orientation  of  Rx  antenna 

variable  (depending  on  direction  of  strong¬ 
est  incoming  signals) 

Table  9.  Urbana  input  parameters  for  simulation  of  angle  diversity  on  receive 

Figure  49  shows  the  five  waypoints  of  the  flight  path  of  a  transmitting  UAV.  The 
waypoints  are  the  locations  where  transmission  is  simulated.  The  signal  strength  coverage 
plots  ofa/L/2  dipole  transmitting  at  location  Tx  #1  to  Tx  #5  are  shown  in  Appendix  D. 

In  investigating  the  perfonnance  of  angle  diversity,  the  receiver  location  (Rx)  at 
(145  m,  52  m,  1.5  m)  uses  the  directional  antenna  pointed  towards  the  strongest  incoming 
signal  ray  for  each  transmitting  location.  The  location  Tx  #3  was  the  same  transmitter  lo¬ 
cation  as  for  the  other  diversity  techniques,  and  therefore  this  point  was  used  to  compare 
performance  with  the  other  diversity  techniques.  The  data  for  the  rest  of  the  transmitting 
locations  are  shown  in  Appendix  D. 
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Figure  49.  Transmit  locations  for  the  UAV  flying  across  the  city 

Figure  50  shows  the  angle-of-arrival  (AO A)  data  for  a  dipole  receiver  at  location 
(145  m,  52  m,  1.5  m)  with  the  dipole  transmitter  at  Tx  #3.  The  vectors  in  Figure  50(a) 
point  to  the  directions  of  the  arriving  rays  at  the  receiver.  Figure  50(b)  gives  the  power 
level  of  each  arriving  ray,  while  Figures  50(c)  and  50(d)  give  the  angles  ^  and  6  of  each 
ray.  Appendix  B  shows  the  representation  of  angles  </>  and  6  in  the  coordinate  system 
used  in  Urbana.  From  these  four  figures,  it  can  be  deduced  that  the  strongest  signal  is  Ray 
Number  #15,  and  it  has  an  AOA  of  (f)  =  45°  and  6  =  -38°. 
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Figure  50.  Signal  rays  received  by  dipole  receiver  at  location  (145  m,  52  m,  1.5  m) 

when  transmitting  from  Tx  #3 


With  the  AOA  data,  the  directional  antenna  was  set  to  point  towards  the  direction 
of  the  highest  power  ray.  The  simulated  results  are  shown  in  Figure  5 1 .  By  comparing 
with  Figures  50(a)  and  5 1(a),  it  is  observed  that  with  the  directional  antenna,  there  are 
fewer  numbers  of  arriving  rays  at  the  receiver  as  compared  to  the  dipole  receiver.  This  is 
primarily  due  to  the  fact  that  the  antenna  radiation  pattern  in  its  rear  hemisphere  is  zero. 
This  is  good  for  rejecting  interference  coming  from  the  rear  direction,  thus  a  good  signal- 
to-interference  ratio.  Looking  at  Figures  50(b)  and  5 1(b),  it  is  also  observed  that  there  is 
improvement  in  the  received  signal  power  for  the  directional  antenna  receiver. 
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Figure  5 1 .  Signal  rays  received  by  directional  antenna  at  location  (145  m,  52  m,  1 .5 

m)  when  transmitting  from  Tx  #3 


In  a  more  detailed  analysis,  Urbana  aggregates  the  signal  strength  of  each  arriving 
ray  and  computes  the  total  received  power  and  outputs  the  result  to  a  couple  fde.  Table 
10  compares  the  total  received  power  at  location  (145  m,  52  m,  1.5  m)  for  the  dipole  and 
directional  antenna  receiver  for  transmit  positions  Tx  #1  to  Tx  #5.  It  is  noted  that  at  the 
deep  fading  location  of  (145  m,  52  m,  1.5  m),  a  directional  antenna  employing  angle  di¬ 
versity  would  be  able  to  provide  a  significant  improvement  in  the  received  signal  power 
(from  1 1  dB  to  23  dB). 


Transmitter 

Received  Power 

Antenna 

Direction 

{.m 

Received 

Power 

Improvement 

Dipole 

Directional 

Antenna 

Tx  #1 

-63.68  dBm 

-51.65  dBm 

(2°, -55°) 

12.03  dB 

Tx  #2 

-66.88  dBm 

-55.61  dBm 

(17°, -52°) 

11.27  dB 

Tx  #3 

-86.96  dBm 

-63.27  dBm 

(45°, -38°) 

23.69  dB 

Tx  #4 

-69.14  dBm 

-48.91  dBm 

(52°, -53°) 

20.23  dB 

Tx  #5 

-67.01  dBm 

-47.00  dBm 

(68°,  -56°) 

20.01  dB 

Table  10.  Total  received  power  for  a  dipole  and  directional  antenna  receiver 
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D.  SUMMARY  OF  FINDINGS 

This  chapter  presented  various  simulations  implementing  the  techniques  of  (1) 
spatial  diversity,  (2)  polarization  diversity,  and  (3)  angle  diversity,  for  UAV  urban  data 
links. 

Table  1 1  shows  a  comparison  of  the  received  signal  power  at  receive  location 
(145  m,  52  m,  1.5  m)  and  transmitting  location  Tx  #3  for  each  diversity  technique  that 
was  investigated  in  this  thesis.  It  is  obvious  that  angle  diversity  has  the  best  performance 
as  compared  to  the  rest,  while  spatial  diversity  with  two  dipoles  also  has  shown  a  signifi¬ 
cant  increase. 


Receiver  at  deep  fading  location 
(145  m,  52  m,  1.5  m) 

Transmitter  at  location 
Tx  #3 

Received 

Signal 

Strength 

One  vertical  dipole  (No  diversity) 

One  vertical  dipole 

-86.96  dBm 

Two  vertical  dipoles  (Spatial  diversity) 

Two  vertical  dipoles 

-81.19  dBm 

Two  horizontal  dipoles  (Spatial  diversity) 

Two  horizontal  dipoles 

-83.19  dBm 

Dual-polarized  dipoles  (Polarization  diver¬ 
sity) 

Dual-polarized  dipoles 

-87.98  dBm 

Two  vertical  dipoles  (Spatial  diversity) 

Dual-polarized  dipoles 

-82.35  dBm 

Two  horizontal  dipoles  (Spatial  diversity) 

Dual-polarized  dipoles 

-84.44  dBm 

Directional  antenna  (Angle  diversity) 

One  vertical  dipole 

-63.27  dBm 

Table  1 1 .  Comparison  of  perfonnance  for  different  diversity  techniques 

For  spatial  diversity,  it  was  demonstrated  that  diversity  on  receive  can  improve 
the  received  signal  strength.  For  example,  with  two  vertical  receive  dipoles  there  is  an 
improvement  of  about  6  dB  at  the  location  of  deep  fading  (145  m,  52  m,  1.5  m),  com¬ 
pared  to  the  case  where  no  diversity  is  used.  The  SNR  improvement  is  3  dB  with  EGC  at 
the  receivers.  Spatial  diversity  with  vertical  dipoles  has  shown  better  performance  over 
spatial  diversity  with  horizontal  dipoles.  The  vertical  case  has  a  smaller  number  of  sig¬ 
nificant  signals  fades  as  compared  to  the  horizontal  case.  The  exception  is  if  the  receivers 
are  located  below  the  transmitter.  Then  horizontal  dipoles  would  provide  better  perfonn¬ 
ance.  Further  improvement  may  be  possible  using  more  than  two  antennas  on  receive,  but 
this  comes  at  the  expense  of  increased  size  and  complexity. 
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With  regard  to  space  limitation,  polarization  diversity  is  a  good  option  to  con¬ 
sider.  A  few  combinations  of  polarization  and  spatial  diversity  were  computed  for  com¬ 
parison.  Polarization  diversity  on  transmit  and  spatial  diversity  on  receive  has  shown  the 
best  performance  of  all  the  combinations.  For  example,  with  dual-polarized  dipoles  on 
transmit  and  spatial  diversity  on  receive,  there  is  a  signal  strength  improvements  of  about 
2.5  dB  and  4.5  dB  for  horizontal  and  vertical  dipoles  on  receive,  respectively.  With  EGC 
at  the  receivers,  the  SNR  improvements  are  -0.5  dBand  1.5  dB,  respectively.  This  is  for 
the  location  of  deep  fading  and  relative  to  the  no  diversity  case.  The  overall  perfonnance 
is  also  not  as  good  as  the  spatial  diversity  technique. 

For  angle  diversity,  the  results  have  shown  that  it  has  much  higher  received  signal 
strength  relative  to  both  the  spatial  and  polarization  diversity  techniques.  For  example, 
when  compared  to  spatial  diversity  on  both  receive  and  transmit  (vertical  cases),  at 
transmitting  location  of  Tx  #3,  the  signal  strength  at  the  deep  fading  location  improved 
by  about  23  dB.  When  angle  diversity  is  compared  with  polarization  diversity  on  both  re¬ 
ceive  and  transmit,  the  improved  received  signal  strength  is  about  25  dB.  In  addition,  the 
angle  diversity  provides  for  good  signal-to-interference  ratio  at  the  receiver.  However, 
the  disadvantage  of  angle  diversity  is  that  the  receiver  would  need  to  be  equipped  with 
the  capability  (i.e.,  an  adaptive  antenna)  to  detect  the  direction  of  the  strongest  arriving 
signal  rays  and  then  to  turn  the  antenna  beam  in  that  direction. 
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VI.  CONCLUSIONS  AND  FUTURE  WORK 


A.  CONCLUSIONS 

The  urban  environment  complicates  military  communications.  The  challenge  for 
MOUT  is  to  achieve  effective  tactical  communications  in  locations  where  radiowave 
propagation  is  subjected  to  multipath  and  fading  because  of  buildings  and  other  structures 
obstructing  the  signals.  UAVs  and  diversity  techniques,  however,  can  be  used  to  mitigate 
these  effects  and  improve  the  reliability  of  communications. 

In  this  thesis,  three  diversity  techniques  were  investigated  to  assess  their  perform¬ 
ance  in  an  urban  environment.  They  were  (1)  spatial  diversity,  (2)  polarization  diversity, 
and  (3)  angle  diversity.  These  techniques  were  simulated  for  a  UAV  platfonn  operating 
over  a  small  city  model.  The  findings  can  be  summarized  as  follows.  Angle  diversity  has 
shown  the  most  significant  improvement  in  received  signal  strength  as  compared  to  the 
other  two  techniques,  but  the  implementation  of  such  directional  antennas  would  be  com¬ 
plex.  Spatial  diversity  on  receive  has  shown  better  overall  improvement  as  compared  to 
polarization  diversity.  Further  improvement  is  possible  with  more  spatially  uncorrelated 
antennas,  but  its  implementation  could  be  limited  by  the  availability  of  space.  Polariza¬ 
tion  diversity  on  transmit  and  receive  has  shown  good  received  signal  strength  only  at 
some  locations  over  the  spatial  diversity  technique.  However,  a  combination  of  polariza¬ 
tion  diversity  on  transmit  and  spatial  diversity  on  receive  has  better  overall  performance. 
An  advantage  of  polarization  diversity  is  that  the  two  orthogonal  linear  polarizations  can 
usually  be  obtained  with  a  single  antenna.  This  is  especially  appealing  for  UAV  applica¬ 
tions,  where  space  on  the  platform  is  limited. 

This  research  has  demonstrated  the  potential  of  implementing  diversity  techniques 
for  improving  UAV  data  links  in  MOUT.  It  should  be  noted  that  the  simulation  per¬ 
formed  did  not  cover  all  possible  system  and  environment  parameters.  A  limited  number 
of  transmitter-receiver  geometries  were  examined.  However,  care  was  taken  to  select 
cases  that  included  direct  LOS  (strong  signals)  and  non-LOS  shadows  (faded  signals). 
This  allows  a  comparison  of  various  techniques  under  a  range  of  conditions,  so  that  some 
general  conclusions  could  be  reached. 
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B.  FUTURE  WORK 

1.  Full  Spectrum  Modeling  and  Simulation 

A  detailed  profding  of  an  exact  operational  urban  area,  e.g.,  downtown  Baghdad, 
would  be  useful.  For  example,  details  like  lamp  posts,  trees  and  roads  can  be  modeled 
and  simulated  to  examine  their  impact  on  the  diversity  techniques.  A  range  of  frequencies 
relating  to  different  types  of  UAV  systems  can  also  be  simulated.  This  would  provide  for 
a  full  spectrum  comparison  of  the  diversity  techniques,  to  include  frequency  diversity. 

2.  Verification  and  Validation  of  Models 

With  detailed  profding  and  modeling,  it  would  be  possible  to  compare  the  model 
with  measured  real-world  data.  For  example,  a  900-MHz  VTOL  UAV  system  using  spa¬ 
tial  diversity  on  receive  can  be  prototyped  for  testing  in  an  actual  environment  to  verify 
and  validate  the  results  of  the  simulation.  This  would  provide  a  feedback  on  the  model’s 
accuracy  and  reliability,  and  thus  enhancements  to  the  model,  if  required. 

3.  Linkages  to  Operation  Research  Models 

There  is  a  possibility  of  integrating  the  modeling  and  simulation  results  in  Urbana 
to  the  existing  C4SIR  operation  models  from  Air  Force  Agency  for  Modeling  and  Simu¬ 
lation  (AFAMS),  and  Army  Model  and  Simulation  Office  (AMSO).  The  idea  is  to  ex¬ 
plore  integration  of  the  findings  into  the  existing  operational  models  for  more  insightful 
pursuit. 


4.  Investigation  of  Channel  Coherence  Bandwidth  and  Coherence  Time 

Further  study  of  the  urban  radiowave  propagation  channel  can  be  extended  to 
analyze  its  channel  coherence  bandwidth  and  time  characteristics.  For  example,  Urbana 
can  be  used  to  compute  the  RMS  delay  spread  of  the  received  signals,  which  can  be  used 
to  compute  the  coherence  bandwidth  and  coherence  time  of  the  channel. 

The  coherence  time  information,  which  defines  the  “static-ness”  of  the  channel,  is 
useful  when  considering  the  Doppler  effects  due  to  motions  of  transmitter,  receiver,  and 
surrounding  objects.  Coherence  bandwidth  information  gives  the  maximum  frequency 
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difference  of  the  signals  that  are  still  strongly  correlated  in  amplitude.  This  limits  the  sig¬ 
nal  bandwidth  [15]. 

5  SDMA  and  Smart  Antenna  Modeling 

For  further  simulation  of  angle  diversity  or  SDMA  for  urban  UAV  data  links, 
adaptive  antenna  receivers  can  be  modeled.  An  adaptive  antenna  would  dynamically  ad¬ 
just  its  radiation  pattern  to  minimize  the  effects  of  noise,  interference,  and  mutipath  [18]. 
The  use  of  an  adaptive  array  to  maximize  SNR  would  result  in  suppression  of  noise  jam¬ 
ming  signals.  Using  an  adaptive  algorithm,  e.g.,  the  least-mean-square  (LMS)  algorithm, 
simulations  can  be  run  to  analyze  the  performance. 

6  Simulation  of  Jamming 

The  effectiveness  of  the  diversity  techniques  against  communications  jamming 
can  be  simulated.  For  example,  a  transmitting  jammer  antenna  can  be  simulated  at  a  cer¬ 
tain  location,  and  then  the  received  jamming  signals  at  various  receiving  locations  can  be 
computed.  Together  with  the  no  jamming  received  signals,  the  jamming-to-signal  ratio 
(JSR)  can  be  calculated.  With  the  various  diversity  techniques  simulated  at  the  receivers, 
different  JSR  can  be  computed  and  compared. 
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APPENDIX  A.  URBANA  CODES 


The  following  are  listings  of  typical  Urbana  user  input  file  and  reprocessing  input 
file  of filename,  ur  input  and  filename.  urrp_input,  respectively: 


-  input  Urbana  v  2.5 

# 

$  ******************************** 

#  A - scatterer  file, length  &  freq 

ft  ******************************** 

#  -  name  of  scatterer  file  in  ACAD  format  (e.g.  wall. facet) 

citywgrnd. facet 

#  -  length  unit:l=inch,  2=cm,  3=meter,  4=mm,  5=mil 

3 

#  -  uniform  freq  (GHz) :  start  freq,  end  ,  nstep 

#  (nstep=0  means:  just  do  first  freq.  CAUTION:  antenna  patterns  are 

#  assumed  to  be  indep.  of  freq  and  is  calculated  at  end  freq) 

0.9  0.9  0 

# 

ft  ********************************* 

#  B -  Antenna  Description  and  List 

ft  ********************************* 

# 

#  - Enter  method  of  describing  antennas. 

#  (1  =  here,  2  =  file) : 

1 

#  - If  described  in  file,  enter  file  name: 

city . antenna 

#  If  described  here,  fill  in  sections  Bl,  B2,  B3 . 

#  If  described  in  file,  use  dummy  data  in  sections  Bl,  B2,  B3 

#  (specify  one  dummy  antenna  type,  dummy  antenna  origin, 

#  and  one  dummy  item  in  antenna  list)  . 

# 

ft  ************************ 

#  Bl :  Define  Antenna  Types 
^  ************************ 

# 

#  Two  lines  for  each  type. 

#  Linel:  type  ID,  ant  code 

#  Line2 :  parameters 

# 

#  Type  ID  must  start  from  1  and  increment  by  1  thereafter 

# 

#  Ant  Code  meaning  parameters 

#  -  -  - 

#  1  pattern  file  filename (ascii) 

#  2  dipole  length (real) 

# 
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#  Antenna  Types  list: 

# 

#  Enter  number  of  antenna  types: 

1 

#  Type  #1 
1  2 

0.1667 

# 

$  *********************************************** 

#  B2 :  Enter  origin  of  antenna  coord  in  main  coord 
$  *********************************************** 

# 

0.  0.  0. 

# 

$  *********************** 

#  B3 :  Create  Antenna  List 
$  *********************** 

# 

#  Three  lines  for  each  antenna. 

#  Linel:  Type  ID,  location  (x,y, z),  power  (watts),  phase (deg) 

#  Line2 :  Local  x-axis  in  main  coord. 

#  Line3:  Local  z-axis  in  main  coord. 

# 

#  Enter  number  of  antennas: 

1 

# 

#  Antenna  #1 

1  -10  -55  150  1.  0. 

1.  0.  0. 

0.  0.  1. 

# 

$  ********************** 

#  C - Observation  points 

ft  ********************** 

#  Observation  points  defined  with  respect  to  main  coord,  system  7. 

#  Enter  method  of  specifying  list  of  points. 

#  (1  =  here,  2  =  file) : 

2 

#  -  If  points  are  listed  here,  enter  number  of  points  (kobtot) : 

1 

#  If  listed  here  (1  above).  List  xyz  of  points  in  main  coord  7 

#  (one  point  at  a  line) .  If  2  above,  include  one  dummy  line. 

1.  2.  -11.00 

#  -  If  points  listed  in  file  (2  above),  enter  name  of  file. 

observl5 . obv 

#  Include  direct  Tx  to  observer  contribution. 

#  If  you  turn  on  the  direct  contribution  from  the  transmitter  to  the 

#  observation  point,  computed  result  will  be  the  total  field,  which 
is 

#  the  incident  +  scattered  field.  For  propagation  analysis,  this 
is 

#  the  preferred  setting.  Otherwise,  the  result  only  includes  the 

#  scattered  field. 

# 
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#  Include  direct  contribution  from  transmitter  to  observation  point 
(rx) 

#  (1  =  yes,  0,2  =  no) : 

1 

#  -  Compute  received  power  into  Rx  antenna. 

#  Urbana  always  computes  field  levels  at  the  observation  point. 

#  If  you  specify  an  Rx  antenna,  Urbana  will  also  compute  the  re¬ 
ceived 

#  power  and  record  the  results  in  the  (runname) .couple  file. 

#  This  causes  a  moderate  but  slow-down  when  using  the  SBR  method 
(below) . 

# 

#  Include  Rx  antenna  (1  =  yes,  0,2  =  no)  : 

1 

#  -  Rx  antenna  specification 

#  Remaining  entries  in  Section  C  can  be  ignored  if  not  including 

#  an  Rx  antenna. 

#  Enter  antenna  type  (1  =  pattern  file,  2  =  dipole) : 

2 

#  Each  antenna  type  requires  additional  parameters. 

#  List  of  expected  parameters  follows.  Choose  one. 

# 

#  Type  Description  Expected  Parameter (s) 

#  1  Pattern  File  File  Name  (e.g.,  beam.antpat) 

#  2  Dipole  Length  (in  prevailing  unit) 

# 

#  Enter  parameter (s)  on  next  line: 

0.1667 

#  -  Rx  antenna  orientation 

#  Enter  local  x-axis  of  Rx  in  global  coordinates 

1.  0.  0. 

#  Enter  local  z-axis  of  Rx  in  global  coordinates 

0.  0.  1. 

# 

$  ***************************** 

#  D - Theoretical  consideration 

$  ***************************** 

#  -  Choose  method  of  computation 

#  0  =  compute  fields  in  the  ABSENCE  of  the  scatterer 

#  1  =  compute  fields  by  SBR 

#  2  =  compute  fields  by  GO 

2 

#  -  If  SBR,  select  a  PO  integration  scheme  at  bounce  points 

#  1  =  do  integration  at  first  &  last  bounce  points  only 

#  2  =  do  so  at  all  bounce  points  (GTD  formulation) 

1 

#  -  Edge  diffraction 

#  SBR  can  be  enhanced  with  PTD  edge  diffraction. 

#  GO  can  be  enhanced  with  GTD  edge  diffraction. 

#  Add  edge  diffraction  (0,2=no,  1=ILDC  (SBR  or  GO),  3=UTD  (GO  only) 
3 

#  -  If  edge  diffraction  switched  on,  enter  name  of  edge  file 

#  (e.g.,  wall. edge  or  dummy  if  edge  not  included) . 

city . edge 
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#  -  Choose  method  of  ray  launch 

#  1  =  by  (baby)  facet,  achieving  a  uniform  first  bounce  surface  den¬ 
sity 

#  2  =  uniform  angular  distribution  (burst  launch) 

#  (If  computation  by  GO,  must  select  2  =  burst  launch) 

2 

#  -  If  ray  launch  by  (baby)  facet  (1  above),  enter  ray  density: 

#  #  rays/wavelength  (normally  5-10) 

5. 

#  -  If  burst  ray  launch  (2  above) ,  enter  angular  interval  (deg) . 

#  (Typically  0.25  -  2.0  deg) 

2. 

#  -  max  permissible  ray  bounces  (normally  5-10) 

7 

#  -  max-voxdepth  =  max  depth  of  BSP  tree  (normally  20) 

#  max-voxl  =  max  facets  in  each  voxel (normally  10) 

#  (Larger  voxdepth  &  smaller  voxl  lead  to  faster  ray-tracing 

#  but  more  computer  memory) 

20,  10 

#  -  ICOAT  for  absorbing  facets 

888 

#  -  IQMATRIX  for  divergence  factor 

#  1  =  calculated  by  Q-matrix 

#  2  =  ignored  except  for  the  spherical  wave  spread 

2 

#  -  IF  using  Q-matrix,  name  target  curvature  file(e.g.  wall.curv) 

dummy . curv 

#  -  IPEC=1  if  all  pec,  =2  if  coating  present 

2 

#  -  For  PEC  scatterer,  give  the  magnitude  of  reflection  coeff 

#  (use  1.0  for  ideal  PEC,  use  less  for  rough  PEC--fudging) 

1.0 

#  -  IF  PEC,  the  rest  coating  info  is  dummmy 

#  -  material  reflection  is  done  through  a  look-up  table 

#  specify  the  freq  interval  in  GHz  for  the  table  e.g.  0.25 

#  (dummy  if  input  freq  less  than  51) 

0.2 

'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

E - coating  material 

'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

-  number  of  materials 

(NOT  including  pec,  which  is  identified  by  ICOAT=0) 

(NOT  including  absorbing  facets:  ICOAT=28  or  888) 

(If  3  material,  urbana  reads  only  ICOAT=l-3) 

3  < - NCOTOT 

-  for  each  material,  identify  its  boundary  type: 

iboundary  =  1  if  impedance  boundary 

2  if  layered  slabs  with  PEC  backing 

3  if  penetrable  layered  slabs  with  air  backing 

4  if  penetrable  FSS  Tx/Refl  table  supplied  by  user 

5  if  same  as  2  except  using  freq-dep  ramlib.d 

6  if  antenna  refl  table  supplied  by  user 

7  if  layers  over  semi-infinite  ground 
for  each  material,  given  info  by  following  the  templates 
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IC0AT=1  A 
iboundary 


3 

-  number  of  layers  over  air  backing 

(1st  layer  is  farthst  fr  incid  field  and  innermost) 

1 


-  thick, epsilon (c) , mu (c) , resistivity (ohm) 

0.30000  (10.1,0.5)  (1.0, 0.0)  l.e+30 

IC0AT=2 


-  iboundary 

2 


-  number  of  layers  over  PEC  backing 

(1st  layer  is  farthst  fr  incid  field  and  innermost) 

2 


-  thick, epsilon (c), mu (c) , resistivity (ohm) 

0.0300  (2.5,-0.000)  (1.6,-0.000)  l.e+30 

0.0500  (3.0,-0.000)  (1.0,-0.000)  l.e+30 

ICOAT=3 


-  iboundary 

7 


-  number  of  layers  over  half-space  (semi-infinite  ground) 

(1st  layer  is  farthst  fr  incid  field  and  innermost) 

1 


- thick,  epsilon (c), mu (c) , resistivity (ohm) 

0.1000  (3.000,-0.0)  (1.0, -0.0)  l.e+30 

-  epsilon (c) , mu (c)  of  semi-infinite  ground 

(3.00,-0.00)  (1.0, -0.00) 


(End  of  regular  input  file.  Leave  a  few  blank  lines) 


'OPTIONAL  ADVANCE  FEATURES'  (Do  not  change  letters  in  quotations) 

#  The  line  above  must  be  placed  at  the  end  of  the  regular  urbana 

#  input.  Advance  features  are  designed  for  special  applications  or 

#  for  testing  codes.  They  are  not  needed  by  general  usages. 

#  - 

#  ADVANCE 1 :  ADD  GTD-TYPE  BLOCKAGE  CHECK 

#  - 

#  In  regular  urbana  computation,  blockage  check  is  mostly  done  by 

#  PTD  principle.  For  interior  scattering  in  a  confined  region,  use  of 

#  GTD  principle  may  be  more  appropriate. 

#  Option  to  use  GTD  principle:  l=yes,  2=no  (regular  case) 

2 

#  - 

#  ADVANCE 2 :  SIMPLE  TERRAIN  BLOCKAGE  MODEL 

#  - 

#  For  GO  method,  terrain  generates  100%  blockage,  and  blocked  rays 
leave 

#  no  energy  behind  a  hill.  With  this  feature,  LOS  rays  and  UTD  edge 

#  diffraction  rays  can  pass  through  terrain,  with  some  attenuation. 

#  Attenuation  is  measured  in  dB  per  hill.  Each  hill  is  identified 

#  by  two  passages  through  two  terrain  facets. 

#  Can  only  be  used  with  GO  method  (and  UTD  edge  option) . 

#  Use  simple  terrain  model:  1  =  yes,  2  =  no  (regular  case) 

2 

#  Enter  coating  code  range  of  terrain  facets  (e.g.,  1,  2): 

1  1 
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#  Enter  amount  of  attenuation  per  hill  (dB,  >  0) : 

5. 

#  - 

#  ADVANCE 3 :  APPROXIMATE  DOUBLE  DIFFRACTION  MODEL 

#  - 

#  For  GO  +  UTD  method,  only  single  diffraction  is  considered. 

#  With  this  feature,  double  diffraction  is  approximated  by  identifying 

#  surfaces  which  block  the  single  diffraction,  such  as  building  walls. 

#  If  one  or  two  facets  block  the  path  from  the  single  diffraction  point 

#  to  the  transmitter,  the  diffraction  is  still  included,  but  with  at¬ 
tenuation  . 

#  Works  best  if  "diffracting  facets",  marked  by  their  coating  code,  are 

#  always  associated  with  enclosed  structures  with  well  defined  edges. 

#  Use  double  diffraction  model:  1  =  yes,  2  =  no  (regular  case) 

2 

#  Encounter  coating  code  range  of  diffracting  facets  (e.g.,  5,  10) : 

2  2 

#  Enter  amount  of  attenuation  for  second  diffraction  (dB,  >  0); 

10. 

#  - 

#  ADVANCE 4 :  ACCELERATION 

#  - 

#  For  large  scenes,  run  time  grows  both  with  the  number  of  field 

#  observation  points  and  the  number  of  edges.  Normally,  all  combina¬ 
tions 

#  of  lit  edges  and  observation  points  are  considered.  This  feature 

#  accelerates  the  processing  by  limiting  the  scope  of  considered  edge 

#  interactions  to  region  around  the  LOS  path  from  the  transmitter 

#  to  the  observation  point.  For  example,  to  run  a  5  km  by  5  km  scene, 

#  one  may  choose  a  250  m  interaction  radius.  For  each  observation 

#  point,  edges  are  ignored  that  lie  outside  an  ellipse  whose  foci  are 
the 

#  Tx  and  the  observation  point  and  whose  major  axis  is  the  LOS  distance 

#  plus  500  m  (radius  x  2) . 

#  This  feature  can  also  be  used  to  automatically  filter  edge  files 

#  whose  domain  far  exceeds  the  domain  of  observation  points. 

#  Only  use  this  feature  for  terrestrial  simulations  where  the  scene 

#  is  nominally  parallel  to  the  x-y  plane. 

# 

#  Use  large  scene  acceleration:  1  =  yes,  2  =  no  (regular  case) 

2 

#  Enter  radius  of  interaction 
250. 

#  - 

#  ADVANCE 5 :  MULTI-DIFFRACTION 

#  - 

#  Substitute  for  Adv.  #3.  Uses  ray  rubber-banding  algorithm  to  find 

#  path  from  transmitter  to  receiver. 

#  Can  only  be  used  with  GO.  Cannot  be  used  in  conjunction  with  Adv.  #3. 

#  If  UTD  switched  on  above,  will  take  measures  not  to  double  count 

#  single  diffraction  mechanisms. 

#  Use  multi-diffraction  model:  1  =  yes,  0,2  =  no 
2 

#  Enter  coating  code  range  of  diffracting  facets  (e.g.  5,  10) : 

2  2 
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#  Enter  maximum  number  of  rubber-band  points  (  also  used  in  Advance6  ) 
1 

#  Check  multiple  crawl  planes  instead  of  just  vertical  one:  1  =  yes, 
0,2  =  no 

0 

#  - 

#  ADVANCE 6 :  REFLECTION-DIFFRACTION 

#  - 

#  If  UTD  switched  on  above,  will  take  measures  not  to  double  count 

#  single  diffraction  mechanisms. 

#  Use  reflection-diffraction  model:  1  =  yes,  0,2  =  no 
2 

#  Do  more  than  just  single  diffractions:  1  =  yes,  0,2  =  no 

#  Allow  rubber-banding  to  both  transmitter  and  receiver:  1  =  yes,  0,2 
no 

1  0 

#  Choose  crawl  plane  selection  mode:  0  =  always  vertical,  1  =  initial 
edge, 

#  2  =  adaptive  from  edge  to  edge 
1 

#  - 

#  ADVANCE 7 :  GREEN'S  FUNCTION  (GF)  FILE 

#  - 

#  By  default,  for  SBR  and  no-target  methods,  a  GF  file  IS  NOT  produced 

#  Also,  by  default,  for  GO,  a  GF  file  IS  produced. 

#  Use  this  feature  to  explicitily  activate  or  de-activate  generation 

#  of  the  GF  file,  which  is  needed  by  the  re-processor  for  its  activi¬ 
ties  . 

#  Activate  GF  file:  0  =  no,  1  =  yes,  2  =  default  activation  behavior 
0 

#  If  yes,  enter  buffer  scale  factor.  Increasing  scale  factor  reduces 

#  the  number  of  GF  file  dumps  to  disk  during  a  run,  but  costs  memory. 

#  Recommend  2-5  for  GO  method,  1  for  no-target  method, 

#  and  100  -  10000  for  SBR  method. 

2 


-  input  Urbana  rp  v  1.2 

# 

$  ******************************** 

#  A - scatterer  file, length  &  freq 

ft  ******************************** 

#  -  name  of  scatterer  file  in  ACAD  format  (e.g.  wall. facet) 

citywgrnd. facet 

#  -  length  unit:l=inch,  2=cm,  3=meter,  4=mm,  5=mil 

3 

#  -  uniform  freq  (GHz) :  start  freq,  end  ,  nstep 

#  (nstep=0  means:  just  do  first  freq.  CAUTION:  antenna  patterns  are 

#  assumed  to  be  indep.  of  freq  and  is  calculated  at  end  freq) 

0.9  0.9  0 

# 

ft  ********************************* 

#  B -  Antenna  Description  and  List 

$  ********************************* 

# 
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#  - Enter  method  of  describing  antennas. 

#  (1  =  here,  2  =  file) : 

1 

#  - If  described  in  file,  enter  file  name: 

cityant . antenna 

#  If  described  here,  fill  in  sections  Bl,  B2,  B3 . 

#  If  described  in  file,  use  dummy  data  in  sections  Bl,  B2,  B3 

#  (specify  one  dummy  antenna  type,  dummy  antenna  origin, 

#  and  one  dummy  item  in  antenna  list) . 

# 

ft  *******  * ***  *  *******  *  ***  * 

#  Bl :  Define  Antenna  Types 
ft  ************************ 

# 

#  Two  lines  for  each  type. 

#  Linel:  type  ID,  ant  code 

#  Line2 :  parameters 

# 

#  Type  ID  must  start  from  1  and  increment  by  1  thereafter 

# 

#  Ant  Code  meaning  parameters 

#  -  -  - 

#  1  pattern  file  filename (ascii) 

#  2  dipole  length (real) 

# 

#  Antenna  Types  list: 

# 

#  Enter  number  of  antenna  types: 

1 

#  Type  #1 
1  2 

0.1667 

# 

ft  *********************************************** 

#  B2 :  Enter  origin  of  antenna  coord  in  main  coord 

ft  *********************************************** 

# 

0.  0.  0. 

# 

ft  *********************** 

#  B3 :  Create  Antenna  List 

ft  *********************** 

# 

#  Three  lines  for  each  antenna. 

#  Linel:  Type  ID,  location  (x,y,z),  power  (watts),  phase (deg) 

#  Line2 :  Local  x-axis  in  main  coord. 

#  Line3:  Local  z-axis  in  main  coord. 

# 

#  Enter  number  of  antennas: 

1 

# 

#  Antenna  #1 

1  -10  -55  150  1.  0. 

1.  0.  0. 

0.  0.  1. 

# 
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$  *******  *  ***  *  ****  *********  * 

#  C - Greens'  function  files 

$  ************************** 

1 

omnirxdiv . gf 

# 

$  ********************** 

#  D - Reprocessor  Output 

ft  ********************** 

# 

#  Generate  field  data  (1  =  yes,  0  =  no)? 

1 

#  Generate  angle  of  arrival  data  (1  =  yes,  0  =  no)? 

1 

#  -  Compute  received  power  into  Rx  antenna. 

#  Urbana_rp  will  record  the  results  in  the  (runname) .couple  file. 

#  If  angle  of  arrival  data  are  also  being  generated,  will  also 

#  record  the  received  signal  per  contribution  in  that  file. 

# 

#  Include  Rx  antenna  (1  =  yes,  0,2  =  no) : 

1 

#  -  Rx  antenna  specification 

#  Remaining  entries  in  this  section  can  be  ignored  if  not  including 

#  an  Rx  antenna. 

#  Enter  antenna  type  (1  =  pattern  file,  2  =  dipole) : 

2 

#  Each  antenna  type  requires  additional  parameters. 

#  List  of  expected  parameters  follows.  Choose  one. 

# 

#  Type  Description  Expected  Parameter (s) 

#  1  Pattern  File  File  Name  (e.g.,  beam.antpat) 

#  2  Dipole  Length  (in  prevailing  unit) 

# 

#  Enter  parameter (s)  on  next  line: 

0.1667 

#  -  Rx  antenna  orientation 

#  Enter  local  x-axis  of  Rx  in  global  coordinates 

1.  0.  0. 

#  Enter  local  z-axis  of  Rx  in  global  coordinates 

0.  0.  1. 
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APPENDIX  B.  COORDINATE  SYSTEM 


The  coordinates  system  used  in  Urbana  is  shown  in  Figure  52. 


z 


Figure  52.  Urbana  coordinates  system 
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APPENDIX  C.  MATLAB  CODES 


The  Matlab  codes  used  to  generate  the  observation  points  and  antenna  pattern  file 
are  shown  below.  Also,  included  are  the  codes  for  computing  the  received  power  for  two 
dipoles  diversity  on  transmit  with  two  dipoles  diversity  on  receive,  and  no  diversity  on 
receive.  Similar  programs  were  written  for  other  combinations  of  diversity,  but  are  not 
included  here.  A  sample  program  for  computing  the  AOA  data  at  the  fixed  receiving  lo¬ 
cation  for  one  of  the  transmitting  location  are  also  shown. 


A.  Generate  Observation  Points 

%  Observation  Points  Generation 

clc; 

clear; 

step  =  2;  %footprint  of  2  m 
i=  1; 

z  =  1.5;  %height  of  1.5  m 

for  x  =  -300:step:300 
for  y  =  -200:step:200 
M(i,:)=[x,y,z]; 
i=i+l; 
end 
end 


lines=i- 1 

save  observ.obv  M  -ASCII; 


B.  Generate  Antenna  Pattern  File 

%  Antenna  Pattern  File  Generation 

clc; 

clear; 

lamda  =  (3*10A8)/(900*10A6); 

HPBW  =  10;  %  in  degrees 

a  =  58.4  *  lamda  /  (2*HPBW)  %  Uniform  Circular  Aperature  Size  (radius) 
beta  =  2*pi/lamda; 

E0  =  1; 

c  =  j*beta*E0*pi*aA2/(2*pi);  %Contant  term  in  E  field 
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Directivity  =  (4*pi*pi*aA2)/lamdaA2 
DirectivitydB  =  I  O*logl  O(Directivity) 

stepl  =  90;  %Vertical  Steps 
step2  =  180;  %Horizontal  Steps 
d_theta  =  1 8  0/(step  1 ) ; 
d_phi  =  360/(step2); 
k  =  0; 
ip  =  0; 
iq  =  0; 

for  phi  =  0:d_phi:360 

for  theta  =  -90:d_theta:90  %plots  out  the  zero  point 
k  =  k+1; 

arg=beta*a*sin(theta*pi/180); 

f  =  l; 

if  abs(arg)>=le-5 
f  =  2*besselj(l,arg)/arg; 
end 

e_theta  =  cos(phi*pi/180)*c*f; 
e_phi  =  -sin(phi*pi/180)*cos(theta*pi/180)*c*f; 
if  theta  >  90 
etheta  =  0; 
e_phi  =  0; 
end 

if  phi  ==  0 
ip  =  ip+  1; 

etO(ip)  =  20*logl0(abs(e_theta)); 
epO(ip)  =  20*logl0(abs(e_phi)); 
eth(ip)  =  theta; 
end 

if  phi  ==  90 

iq  =  iq+ 1; 

et90(iq)  =  20*logl0(abs(e_theta)); 
ep90(iq)  =  20*logl0(abs(e_phi)); 
eth90(iq)  =  theta; 
end 

A(k,l:4)  =  [real(e_theta),  imag(e_theta),  real(e_phi),  imag(e_phi)]; 
end 
end 

save  antennapat  -ASCII; 
figure(l) 

plot(eth,et0,eth,ep0) 

title('Directional  Antenna  Pattern  (HPBW=10  degrees)') 
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xlabel(’Angle  [degree]’) 
ylabel('Relative  Power  [dB]’) 
figure(2) 

plot(eth90,et90,eth90,ep90) 
xlabel(’Angle  [degree]’) 
ylabel('Relative  Pattern,  dB') 


C.  Computation  of  Received  Power  Data 

%Received  Power  for  Diversity  on  transmit  and  receive 

clc; 

clear; 

Txantdiv=2; 

TF=0; 

i~\d=i'opcn(' city2antx.couple\'r');  %  2  dipoles  transmit  diversity  couple  fde 
lineO=fgetl(fid);  %  to  rid  of  the  header 
while  TF==0 
linel=fgetl(fid); 

obv=sscanf(linel,'%i  %*i  %*f  %*f  %*f); 
antid=sscanf(linel,’%*i  %i  %*f  %*f  %*f);  %  antenna  id 
powerdBm(obv,antid)=sscanf(linel,'%*i  %*i  %*f  %f  %*f);  %  power  value 
TF=feof(fid); 
end 

fclose(fid); 

for  n=  1  :obv- 1  %coherently  adding  at  two  observation  points 

powerdiv  1  (n)=(  1 0A((powerdBm(n,  1 ))/ 1 0))+(  1 0A((powerdBm(n+ 1 , 1  ))/l  0)); 
powerdiv2(n)=(  1 0A((powerdBm(n,2))/ 1 0))+(  1 0A((powerdBm(n+ 1 ,2))/ 1 0)); 
powerdivlBm(n)=10*logl0(powerdivl(n));  %  power  from  Txl  to  2  Rx  div 
powerdiv2Bm(n)=10*logl0(powerdiv2(n));  %  power  from  Tx2  to  2  Rx  div 
powerdivtot(n)=powerdivl(n)+powerdiv2(n);  %  coherent  addition  of  2  Rx 
powerdivtotdBm(n)=  1 0*log  1 0(powerdivtot(n)); 
powertot(n)=(  1 0A((powerdBm(n,  1))/ 1 0))+(  1 0A((powerdBm(n,2))/ 10)); 
powertotdBm(n)=10*logl0(powertot(n));  %  power  from  Txl  &  Tx2  to  1  Rx  (no  div) 
end 


n=l:obv-l; 
hold  on 

plot(powertotdBm(n),’b')  %  2  Tx  and  no  Rx  div 
plot(powerdivtotdBm(n),’r’)  %  2  Tx  and  Rx  div 
xlabel(’Observation  points') 
ylabel('P[dBm]') 
hold  off 
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D.  Computation  of  AOA  Data 

%  AOA  Data  Generation  for  one  receiving  location 

clc; 

clear; 

n=0; 

cnt=0; 

TF=0; 

fi  d=  fo  p  c  n  ('  om  ni9 00M_  Txl  ,aor\'v'); 

while  TF==0 
lineO=fgetl(fid); 

f=sscanf(lineO,'%*s  %*s  %i%*s  %*s  %*s  %i'); 


for  n  =  1  :f(2)  %no  of  rays 
linel=fgetl(fid); 
line2=fgetl(fid); 
line3=fgetl(fid); 

antid=sscanf(linel,’%i  %*i  %*i  %*f  %*f  %*f  %*f  %*f); 

r=sscanf(line3,'(%f,%*f)'); 

im=sscanf(line3,'(%*f,%f)'); 

raypoweri=(rA2  +  imA2);  %for  testing  of  valid  ray  only 
if  (raypoweri  >  0) 
cnt=cnt+l; 

type(f(l),antid,cnt)=sscanf(linel,'%*i  %i  %*i  %*f  %*f  %*f  %*f  %*f);  %ray  type 

bounce(f(l),antid,cnt)=sscanf(linel,'%*i  %*i  %i  %*f  %*f  %*f  %*f  %*f);  %no. 
of  bounces  of  ray 

x(f(l),antid,cnt)=sscanf(linel,'%*i  %*i  %*i  %*f  %*f  %f); 
y(f(l),antid,cnt)=sscanf(linel,'%*i  %*i  %*i  %*f  %*f  %*f  %f); 
z(f(l),antid,cnt)=sscanf(linel,'%*i  %*i  %*i  %*f  %*f  %*f  %*f  %f); 

ph=atan(y(f(  l),antid,cnt)/x(f(  l),antid,cnt)); 

phi(f(l),antid,cnt)=rad2deg(ph); 

%phi_angle  for  one  observation  point,  tx  antenna,  ray 

theta(f(l),antid,cnt)=rad2deg(atan(y(f(l),antid,cnt)/(z(f(l),antid,cnt)*sin(ph)))); 

%  theta  angle  for  one  observation  point,  tx  antenna,  ray 

raypowerf(f(l),antid,cnt)=raypoweri;  %actual  ray  power 
end 
end 

ray(f(l))=cnt;  %no  of  actual  rays  for  the  fixed  receiving  point 
cnt=0; 
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TF=feof(fid); 

end 

fclose(fid); 

figure(l) 

q=l  :ray(f(l));  %get  the  exact  no  of  rays  for  the  fixed  receiving  point 

p=phi(f(l),l,q); 

clear  s; 

s(q)=p; 

plot(q,s,':s') 

title(’AOA  from  Tx#') 

xlabel('Ray  Number') 

ylabel('\phi  Degree') 

figure(2) 
t=theta(f(l),l,q); 
clear  u; 
u(q)=t; 
plot(q,u,':s') 
title(’AOA  from  Tx#') 
xlabel('Ray  #') 
ylabel('\theta  Degree') 

figure(3) 

rp=  1 0*log  1 0((raypowerf(f(  1 ),  1  ,q))*  1 000); 

clear  w; 

w(q)=rp; 

plot(q,w,':s') 

title('Power  of  rays  from  Tx#') 
xlabel('Ray  #') 
ylabel('P[dBm]') 

figure(5) 

gx=gradient(x(f(l),  1  ,q)); 
gy=gradient(y(f(  1),  1  ,q)); 
gz=gradient(z(f(  1),  1  ,q)); 

quiver3(gx,gy,gz,x(f(  1),  1  ,q),y(f(  1),  1  ,q),z(f(  1),  1  ,q),  1 00)  %plot  the  aor  of  rays 
view([70  18]) 

title('Direction  of  rays  from  Tx#') 

xlabel('x') 

ylabel(’y') 

zlabel('z') 
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APPENDIX  D.  SIGNAL  COVERAGE  AND  RAYS  PLOTS  FOR 

ANGLE  DIVERSITY 


The  signal  coverage  and  rays  plots  for  Tx  #1  to  Tx  #5  are  shown  below. 


A.  Signal  Coverage  Plots  for  Tx  #1  to  Tx  #5 
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Figure  54.  Signal  coverage  for  Tx  #2 


Figure  55.  Signal  coverage  for  Tx  #3 
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Figure  56.  Signal  coverage  for  Tx  #4 
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Figure  57.  Signal  coverage  for  Tx  #5 
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B 


Signal  Rays  Plots  for  Tx  #1 
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Figure  58.  Signal  rays  received  by  dipole  receiver  at  location  (145  m,  52  m,  1.5  m) 

when  transmitting  from  Tx  #1 


(a)  Ray  AO  As  plot  (b)  Ray  power  levels 

Figure  59.  Signal  rays  received  by  directional  antenna  at  location  (145  m,  52  m,  1.5 

m)  when  transmitting  from  Tx  #1 
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C.  Signal  Rays  Plots  for  Tx  #2 
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Figure  60.  Signal  rays  received  by  dipole  receiver  at  location  (145  m,  52  m,  1.5  m) 

when  transmitting  from  Tx  #2 
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Figure  61.  Signal  rays  received  by  directional  antenna  at  location  (145  m,  52  m,  1.5 

m)  when  transmitting  from  Tx  #2 
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D.  Signal  Rays  Plots  for  Tx  #4 
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Figure  62.  Signal  rays  received  by  dipole  receiver  at  location  (145  m,  52  m,  1.5  m) 

when  transmitting  from  Tx  #4 
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Figure  63.  Signal  rays  received  by  directional  antenna  at  location  (145  m,  52  m,  1.5 

m)  when  transmitting  from  Tx  #4 
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E.  Signal  Rays  Plots  for  Tx  #5 
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Figure  64.  Signal  rays  received  by  dipole  receiver  at  location  (145  m,  52  m,  1.5  m) 

when  transmitting  from  Tx  #5 
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Figure  65.  Signal  rays  received  by  directional  antenna  at  location  (145  m,  52  m,  1.5 

m)  when  transmitting  from  Tx  #5 
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